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In a class of shape realization processes, the product shape is realized by making the 
original material in the form of a deformable sheet to conform to a relatively rigid 
surface of the tool. Such shape conformance involves an elastic/plastic deformation of 
the material by one or more of the modes of deformation, namely, bending, stretching 
and in-plane shear. As the shape realization here is due to close conformance of one 
surface to another, these processes can be called as shape conforming processes. Press 
brake bending of sheet metal componeirts, stretch forming, deep drawing, spinning, 
composite layup, compression molding with sheet molding compound (SMC), some of 
the preforming techniques used in resin transfer molding and thermoforming are some 
of the examples in this category. 

Tooling is an important aspect in all shape conforming processes. Tooling in the form 
of molds and dies provide a single or multiple base surfaces required for shaping the 
metallic sheet or the combination of reinforcing and matrix materials. If we neglect the 
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the shape variations due to springback and thermal effects, the shape of the product 
obtained in these processes is completely defined by the shape of the mold or die used. 
In an ideal situation the problem of shape realization in these processes reduces to that 
of realization of a tool surface. 

All shape conforming processes suffer from a limitation in the sense that a tool is 
required as a precursor for manufacturing an object. For trial production runs and 
prototyping work, when the shape changes are frequent, a considerable time and money 
is spent- on the development of a tooling. For example, in the case of composites the 
mold design and fabrication often involves a slow and labour intensive stage of master 
preparation. In a typical prototyping environment, w’here the accuracy requirements of 
the product shapes are not very critical, the hard tooling is generally expensive and can 
be substituted by a flexible tooling for a class of component shapes. One such tooling 
which has been proposed for shape conforming processes is discrete surface tooling. 

The discrete surface tooling is based on the concept of discrete approximation to a 
continuous surface of a mold or a die. It consists of a number of closely spaced multiple 
rigid surface tool elements, each of which is a surface element of an expected contour. 
The heights of the tool elements can be adjusted to approximate the desired surface 
shapes. The positioning of each surface tool element can be carried out either manually 
or using a computer control. The discrete nature of the tool can be overcome by 
placing a deformable elastomeric sheet over the surface tool elements. The flexible 
sheet draped over the discrete surface provides with a continuous surface required for 
the tooling applications. The advantage of such a tooling is that it is adjustable. A 
variety of surface shapes can be realized by properly adjusting the heights of surface 
tool elements. 

The concept of discrete surface tooling was first given by Nakajima (1969). The concept 
was used to fabricate complicated dies for press work and electrodes for electrochemical 
machining. Here, one set of surface tool elements form one half of the mold surface 
which is held in a press and the other set of tool elements form a second half of the mold 
whose motion is controlled by the a ram. An extensive research work on the concept 
of discrete surface tooling has been carried out by Hardt et. al at the Laboratory 
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for Maiiufacturing and Productivity at Massachusetts Institute of Technology. They 
propose the use of a variable configuration discrete element die, which is similar to 
Nakajima’s flexible die. Unlike the forming system proposed by Nakajima, they propose 
the use of discrete element die in a closed loop shape control system. The system forms a 
workpiece by an iterative procedure involving forming, measuring and die shape change 
until the shape convergence occurs. 

A modified version of the above scheme of tooling has been proposed as part of the 
present work which is meant for those processes which can provide a single finished 
surface of the product, such as composite layup. The tool in this case consists of a 
rubber-like sheet constrained partially or fully along its edges and lifted by a set of 
surface tool elements deforming the flexible sheet to the required shape. The tooling 
in this case is known as flexible surface tool and the surface tool elements here are 
referred as indentors. One side of the deformed flexible sheet conforms to discrete 
surface formed by the indentors and the other side of the sheet can be used for tooling 
applications. 

The present research is motivated by the need to explore the possibilities of using the 
proposed flexible surface tooling for shape conforming processes in general, and for 
composite layup in particular. Before a flexible surface tooling machine based on the 
proposed concept is built, it is necessary to evaluate such an idea computationally. 
This will help in identifying the the class of shapes which can be realized. The shape 
accuracies which can be obtained with such a tooling is another important aspect. The 
computational simulation of the proposed concept also forms a basis for the successful 
design and the fabrication of the flexible surface tooling machine. 

The successful realization of the desired shapes by a flexible surface tooling depends 
on many tooling parameters. These include the size, the thickness, and the material 
of the flexible elastomeric sheet, the layout of indentors, the contact and the section 
geometries of the indentors etc. The study of these tooling parameters and their effect 
on the shapes obtainable is necessary. An experimental study of these parameters is a 
time consuming and also not cost effective. 
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The overall object.ive of the present work is to carry out the numerical study of the 
flexible surface tooling concept in order to evaluate its applicability to composite layup 
process in a prototyping environment. The specific objective of this research can be 
stated as follow's: 

• To identify the shape conforming processes and to propose the use of discrete 
surface tooling concept to this class of processes in general and composite layup in 
particular. 

A class of manufacturing processes have been identified where the shape realization is 
due to geometric transformation from an original shape to a tool shape. A discrete 
surface tooling which is based on the concept of discrete approximation to a continuous 
surface, has been extended for shape conforming processes in general and composite 
manufacturing processes in particular. 

• A computational study of the realizability of a class of shapes by flexible surface 
tooling. 

A numerical st udy of the realizability of simple and complex shapes by the proposed 
flexible surface tool has been carried out. A computer program which deals with 
the deformation analysis of nonlinear elastic membranes in multiple contact has been 
written and used for this purpose. 

• To study the effect of tooling parameters on the surface shape accuracies. 

The, effect geometric and material properties of the flexible sheet on the shape accura- 
cies has been studied with an objective to identify the range of these parameter for a 
given application. 

• A study of the parameters associated with discrete approximation to a continuous 
surface. 

The tooling parameters related to surface discretization such as the number and the 
layout of indentors, the section and the contact geometries of the indentors play a 
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major role in obtaining accurate surface shapes. The effect of these parameters on the 
accuracy of surface shapes has been studied with an objective to identify the suitable 
range of theses values for fabricating a flexible surface tooling machine. 

• To study the static shape optimization of the flexible sheet in the proposed tooling. 

The static shape of the deformed flexible sheet in the proposed tooling is used as 
a base surface for tooling applications. Apart from the other parameters associated 
with the tooling, the deformed shape of the flexible sheet depends on the heights by 
which indentors are raised. The sensitivity of the static shape to the indentor heights 
is an important aspect. This has been dealt with an objective to improve the shape 
accuracies obtainable. 

• To study the realizability of shapes having negative Gaussian curvature. 

The proposed scheme of tooling has been modified for realizing the shapes having 
negative Gaussian curvature. A numerical study of this scheme of tooling has been 
carried out in order to evaluate the realizability of shapes. 

Chapter 1 of the thesis introduces the shape conforming processes and the importance 
of tooling in these processes. The concept of adjustable discrete tooling is discussed 
and its applicability for shape conforming processes in general is discussed. A modified 
version of the discrete surface tooling namely, flexible surface tooling is proposed for 
shape conforming processes where a single finished surface is desired. The need for a 
computational analysis of this version of tooling is discussed. 

In the proposed tooling system, deformed configuration of the flexible elastomeric sheet 
is used as the mold surface. Hence, the deformation analysis of the flexible sheet forms 
the central idea of the present work. This involves solving a finite elasto-static contact 
problem with multiple contacts. Chapter 2 presents the formulation and the numerical 
scheme used for solving this problem. The chapter also presents the validation of the 
computer code developed for the purpose of analysis. 

Chapter 3 discusses the description of the steps involved in the shape realization proce- 
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dure in terms of a preprocessor, analysis and the post processor. The chapter presents 
the results of the numerical experiments in terms of realization of singly and doubly 
curved surface shapes. The results of the study of various parameters associated with 
the flexible sheet such as geometric and material properties of the sheet, sheet dimen- 
sions etc. on the accuracy of realized shapes are presented. The effect of external 
loading on the realized shapes is also dealt as a part of this chapter. 

Clhapter 4 is concerned with the study of parameters involved in discrete approximation 
of a continuous surface. These include the number of indentors, the indentor layout, 
the section and contact geometries of the indentors etc. The objective of this study 
is twofold. The main purpose of this study is to know the effect of discretization 
parameters in order to improve the accuracy of realizable shapes. The secondary 
objective of such a study is to identify the range of discretization parameters in order 
to fabricate a flexible surface tooling machine. 

Chapter 5 of the thesis deals w'ith the static shape optimization of the flexible sheet 
for the indentor heights. The presence of a flexible sheet on the discrete surface of the 
mold changes the effective mold shape. In such cases, the sensitivity of the static shape 
of the flexible sheet for indentor heights forms an important and interesting aspect of 
the present study. This aspect has been investigated with an objective to minimize the 
error between the desired and realized surface. 

An alternate version of the flexible surface tooling for realizing shapes having both syn- 
clastic and anticlastic regions is discussed. A numerical study of this version of tooling 
forms the contents of Chapter 6. The main objective of this study is to evaluate the 
realizability of shapes having considerable negative Gaussian curvature. A secondary 
objective of this study is to compare this version of tooling with the earlier version 
dealt as a part of Chapter 3. 

The conclusions of this work are presented in Chapter 7. The future scope of the work 
in terms of analysis, design and fabrication of the proposed tooling are discussed in 
detail. The possibility of extending the concepts and results which are the outcome of 
the present work to applications in other areas of engineering is also discussed. 
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Chapter 1 


INTRODUCTION 


1.1 Shape Realization 

In an environment of concurrent engineering, the physical realizability of a designed 
shape by a manufacturing process is of primary concern. The shape realization by 
a manufacturing process often involves achieving the desired shape of a product or a 
component from an original shape by using a proper tool shape and the associated 
machine tool kinematics (Fig. 1.1). Most of the shape realization processes can be 
classified into two major groups, namely, shape generation processes and shape for- 
mation or shape transformation processes (Tavakkoli, 1991; Dhande and Karunakaran, 
1993). 

In a shape generation process, the product shape is realized by the removal of material. 
In general, such material removal requires movement of a cutting tool through the 
volume of a blank (Fig. 1.2). The tool shape as well as the kinematics of machine tool 
play an important role in the process of shape realization. Processes such as milling, 
gear bobbing, electric discharge machining and sheet metal punching can be classified 
as shape generation processes. The realizability of a given shape by shape generation 
processes has been studied extensively by Voruganti (1990), Karunakaran (1993) and 
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Figure 1.1. Shape Realization. 
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Misra (1993). 

A shape forming process, on the other hand, involves a geometric transformation of 
a shape from its original state to a desired state (Fig. 1.2). Such transformations 
often involve a permanent or plastic deformation of the material. In a shape forming 
process, the tool shape plays a major role in realizing the shape of a product. The 
associated machine tool kinematics is relatively simple and its role is to bring the 
original material in contact with the tool surface by deforming it. In a shape forming 
process, the shape transformations can be idealized either at a volume level when they 
are called as bulk forming processes or at a surface level when they are called as sheet 
forming processes. The processes such as forging, rolling and extrusion can be classified 
as bulk forming processes, whereas sheet bending, deep drawing and roll forming are 
some of the processes which can be categorized as sheet forming processes (Alting, 
1982). 


1.2 Shape Conforming Processes 

In a class of sheet forming processes, the shape is realized by making the original 
material in the form of a deformable sheet to conform to a relatively rigid surface 
of the tool. Such shape conformance involves an elastic/plastic deformation of the 
material by one or more of the modes of deformation, namely, bending, stretching and 
in-plane shear. As the shape realization here is due to close conformance of one surface 
to another, these processes can be called as shape conforming processes. In other words, 
in all shape conforming processes the product shape is dictated by the shape of a tool. 
The following paragraphs describe some of the important shape conforming processes 
and the role of tool/mold in shape realization. 

A number of sheet metal forming processes used with metallic materials can be classified 
as shape conforming processes. In these processes, the sheet metal in the form of 
a blank undergoes a permanent deformation by the application of force to conform 
to complicated die shapes. Press brake bending of sheet metal components, stretch 
forming, deep drawing and spinning are some of the examples in this category (Fig. 
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1.3). Neglecting the springback associated with the deformation, the shape of the 
product or a component realized in these processes is completely defined by the shape 
of a tool used. 

A majority of the composite fabrication processes can be classified as shape conform- 
ing processes (ASM Handbook, 1988; Mallick and Newman, 1990; Bader et. al, 1990; 
Potter, 1989). In layup, one of the primary composite fabrication processes, the mul- 
tiple layers of fiber and resin are placed either by hand or a machine into a preshaped 
mold until the desired thickness is reached. This is often followed by the application of 
pressure and heat to cause the part to consolidate and cure. The shape of the product 
obtained here is dictated by the mold shape (Fig. 1.4). Similarly, in one of the versions 
of compression molding process, the mixture of resin and fiber in the form of a sheet of 
ready-to-mold composite material known as sheet molding compound(SMC) is placed 
into one half of the mold cavity mounted in a platen press and is closed with the other 
half. Under the applied pressure and elevated temperature, the material deforms to 
conform to the mold surfaces when it is allowed to cure. The shape of the product ob- 
tained in the process is defined by mold surfaces to which the sheet molding compound 
is made to conform (Fig. 1.5). 

In a resin transfer molding (RTM) or liquid composite molding (LCM) process, the 
shape of the product is defined at the stage of preforming. A preform is an assembly 
of fibrous reinforcement and other inserts that are pre-oriented and formed to near or 
near-net shape serving as a skeleton of the actual structure. In the common preforming 
techniques such as cut-and-place preforming, directed-fiber preforming and stamping 
of thermoformeable materials, the reinforcement material is given a shape of a preform 
using the predefined mold surfaces. In this manner, it is similar to layup process and 
the shape realization is conceptually same as any of the shape conforming processes 
discussed above (Carley et. al, 1989; Steenkamer et. al, 1993). 

Manufacture of fiber reinforced thermoplastic composites using preforms made of com- 
mingled fabrics is another example in the category of shape conforming processes. 
Commingled fabric is composed of yarns containing both the reinforcement and matrix 
in fiber form. A quality of the commingled fabric is its drapability allowing it to con- 
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Figure 1,3. Shape Conforming Processes I. 
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form to mold surfaces having a compound curvature prior to being consolidated into a 
rigid structure (Shukla, 1989; Van West, 1990). The process of draping here involves 
mainly two modes of deformation, namely, bending and in-plane shear. 

Another example of shape conforming process is thermoforming. The process starts 
with heating of a thermoplastic sheet to its softening or sag temperature. The hot 
thermoplastic sheet in rubbery state is placed over a mold and is drawn by vacuum 
or pressure into or over the mold (Fig. 1.6). The material is held in the mold until it 
cools and takes the shape of the mold (Farnham, 1972). 


1.3 Tooling for Shape Conforming Processes 

Tooling is an important aspect in all shape conforming processes. Tooling in the form 
of molds and dies provide a single or multiple base surfaces required for shaping the 
metallic sheet or the combination of reinforcing and matrix materials. In other words, 
tool is a precursor for fabricatioir of every simple or complex shape. If we neglect the 
the shape variations due to springback and thermal effects, the shape of the product 
obtained in these processes is completely defined by the shape of the mold or die used. 
It can be said that, in an ideal situation, the realizability of a shape in all shape 
conforming processes is mainly governed by the tool used and is independent of the 
process. This reduces the problem of shape realization in these processes to that of 
realization of a tool surface. 

The tool design and fabrication for shape conforming processes is a time consuming 
and often a costly affair. This is true even when the number of components required is 
small or the environment is that of prototyping. This is because, many of the processes 
require tools and equipment similar to those used in the actual process to accurately 
simulate the physical properties achievable in the production level components. Typi- 
cally, one spends more time on the design and fabrication of a tooling than that involved 
in the fabrication of component itself. For example, the design and fabrication of a die 
for a free-form sheet metal component is a costly and often experimentally iterative 
process. A number of forming experiments with improved die shapes at each stage are 
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required to converge to a final die shape which can be used for realizing the desired 
component shapes (Hardt and Webb, 1982; Hardt, 1992). The computer aided design 
of dies for such surface shapes calls for the solution of frictional contact problems by 
numerical methods such as FEM in order to compensate for springback. In a proto- 
typing environment, where the number of components required is often limited, such 
numerical simulations or forming experiments are both costly and time consuming. 


1.4 Tooling for Composites 

In the case of composites, the mold design and fabrication is often a slow, labour inten- 
sive and a costly affair. Moreover, the choice of tooling for composites is a complicated 
decision that must take into account many factors. In fact, their design requires the 
consideration of as many factors as are studied in the design of part itself. Presently, 
there are three approaches to the fabrication sequence of the a mold for composites 
(Ham and Molitor, 1989; Strong, 1993). 

In the first approach, which is also a traditional composite tool building process, molds 
are generally made of composite materials such as Glass-Epoxy, Graphite-Epoxy and 
Fiberglass. This approach consists of two stages (Fig. 1.7). First stage involves the 
transfer of tool geometry data to a physical working model. Models are made from 
plaster, Epoxy-faced plaster, Fiberglass-Epoxy, wood and other materials. In a typical 
scheme, full size mylars are used to layout templates. These templates are then used 
to create lofting lines for the splining of either plaster or other materials. In the second 
stage, the master model thus fabricated is used as a base surface for producing a female 
composite mold by wet or prepreg layup process. The availability of room temperature 
curable tooling prepreg systems allow the construction of high quality composite tools 
to be built directly from the master model. 

The above method of master fabrication is a slow and labour intensive process. More- 
over, due to involvement of human element and thermal expansions, the transforma- 
tion of geometric data to a master model is not accurate. In the second approach, 
the masters are produced directly by using common metal working techniques and can 



Chapter 1. Introduction 


12 



Approach I 



Approach II 


Metallic Mold Fabrication 
by 

Machining/Casting/Forming 


Approach III 


Figure 1.7. Fabrication Sequence of a Mold for Composites. 
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subsequently be used for the production of female composite molds (Fig. 1.7). By 
incorporating computer aided design and computer aided manufacturing techniques, it 
is possible to replace an inherently inaccurate hand built model with a master that is 
machined directly from the geometric model. 

In the third approach, the time consuming phase of master preparation is totally elim- 
inated by producing the female molds directly by machining or casting (Fig. 1.7). 
When high volumes or complex multipiece tools are required, molds can be machined 
directly using a CNC machine. In this case, the molds are generally made of Alu- 
minium, Steel or Electroformed Nickel. This shortened tooling sequence provides a 
substantial reduction in lead time. However, for short production runs or prototyping 
work this approach is not economical. Moreover, in many cases it is easier to machine 
the male part of the shape than the female counterpart. 

When the quantity of parts to be produced and the production rate are small, the 
tooling costs would normally be small, unless very high accuracy is required and paid 
for. The machinery for making the parts would also be lower in cost with relatively 
more human craftsmanship involved. For a low volume production, the first approach of 
tool fabrication is still followed in many cases when the molds are made of composites. 
But, this approach of prototyping composite products by preparation of hand-built 
masters, and subsequent composite tooling, followed by layup is a time consuming 
one and hence not productive. There is a need for faster and economical prototyping 
methods to reduce the time involved in marketing a composite product. 


1.5 Flexible Surface Tooling 

All shape conforming processes suffer from a limitation of, requirement of tool as a 
precursor for manufacturing an object. For trial production runs and prototyping 
work, when the shape changes are frequent, a considerable time and money is spent 
on the development of a tooling. Moreover, in a typical prototyping environment, the 
accuracy requirements of the product shapes are not very critical. In such cases, the 
hard tooling is generally expensive and can be substituted by a flexible tooling for 
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a class of component shapes. One such tooling which has been proposed for shape 
conforming processes is discrete surface tooling. 

The discrete surface tooling is based on the concept of discrete approximation to a 
continuous surface of a mold or a die. It consists of, a number of closely spaced multiple 
rigid surface tool elements, each of which is a surface element of an expected contour. 
The heights of the tool elements can be adjusted to approximate the desired surface 
shapes (Fig. 1.8). The positioning of each surface tool element can be carried out either 
manually or using a computer control. The discrete nature of the tool can be overcome 
by placing a deformable elastomeric sheet over the surface tool elements. The flexible 
sheet draped over the discrete surface provides with a continuous surface required for 
the tooling applications. The advantage of such a tooling is that it is adjustable. A 
variety of surface shapes can be realized by properly adjusting the heights of surface 
tool elements. The total time involved in the tool setup is considerably less compared 
to that involved in the development of a hard tool. 

Based on the above concept, a reconfigurable tooling system has been proposed for 
those processes in which two finished surfaces can be obtained such as deep drawing 
and compression molding (Nakajima, 1969; Hardt and Gossard, 1980). Here, one set of 
surface tool elements form one half of the mold surface which is held in a press and the 
other set of tool elements form a second half of the mold whose motion is controlled 
by the a ram (Fig. 1.9). As the surfaces formed by the surface tool elements will 
not be smooth, two flexible sheets of rubber-like material are placed between the the 
component to be formed and each of the mold surfaces respectively. 

A variation of the above scheme of tooling has been proposed here and is meant for 
those processes which can provide a single finished surface of the product, such as 
layup. It consists of a rubber-like sheet constrained partially or fully along its edges 
and lifted by a set of surface tool elements deforming the flexible sheet to the required 
shape (Fig. 1.10). One side of the deformed flexible sheet conforms to discrete surface 
formed by the surface tool elements and the other side of the sheet can be used for 
tooling applications. The tooling in this case is known as flexible surface tool and 
the surface tool elements here are referred as indentors. A relatively small number of 
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Figure 1.8. A Discrete Approximation to a Continuous Surface. 
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indentors are used here and the tooling is meant for those processes where the forming 
forces are considerably less. 

One of the limitations of this scheme of tooling is that component shapes having only 
synclastic regions can be realized. This limitation can be over come by modifying 
the proposed tooling to include the shapes having anticlastic regions. In the modified 
version, the tooling consists of a rubber-like sheet wdth some or all the indentors bonded 
to it. The indentors can be attached to the flexible sheet by mechanical means or, in 
another version, the flexible sheet is held by a set of hollow indentors which act as 
suction cups. However, the nature of tooling restricts the realizable shapes to those 
having small curvature values in this case. 

The two schemes of tooling discussed above namely, the reconfigurable tooling system 
of Nakajima and the flexible surface tooling system proposed here conceptually differ in 
terms of realization of the tool surface. In the first scheme of tooling, the mold surfaces 
are primarily defined by the surface tool elements and the flexible sheet only acts as a 
smoothing material. In this case, the thickness of the sheet is relatively small, of the 
order of 1 to 2 mm and it is subjected to high compressive forces due to sandwiching 
between the component formed and the discrete die surface. On the other hand, in 
the second scheme of tooling the flexible sheet is subjected to tensile forces due to 
indentation by surface tool elements. There is a significant change in the effective die 
shape due to draping of a relatively thick flexible sheet over the discrete surface. In 
this case the mold surface is primarily defined by the deformed configuration of the 
flexible sheet. 


1.6 Motivation of the Present Work 

The present research is motivated by the need to explore the possibilities of using the 
proposed flexible surface tooling for shape conforming processes in general, and for 
composite layup in particular. The successful development of a flexible surface tooling 
will result in considerable reduction of the time involved in prototyping certain class of 
shapes. Moreover, automated nature of the process may relegate the need of some of 
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Figure 1.10. Flexible Surface Tooling. 
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the experts/skilled labours involved in tool design and development. Experiments with 
new shapes and incorporating design changes will be much easier once such a tooling 
is made available. 

Before a flexible surface tooling machine based on the proposed concept is built, it is 
necessary to evaluate such an idea computationally. This will help in identifying the 
the class of shapes which can be realized. The shape accuracies which can be obtained 
with such a tooling is another important aspect. The computational simulation of the 
proposed concept also forms a basis for the successful design and the fabrication of the 
flexible surface tooling machine. 

The successful realization of the desired shapes by a flexible surface tooling depends 
on many tooling parameters. These include the size, the thickness, and the material 
of the flexible elastomeric sheet, the layout of indentors, the contact and the section 
geometries of the indentors etc. The study of these tooling parameters and their effect 
on the shapes obtainable is necessary. An experimental study of these parameters is 
a time consuming and also not cost effective. A theoretical study and computational 
simulation based on finite element analysis is now' possible and is recommended. 


1.7 Overview of Previous Research 

The concept of surface discretization for tooling requirements is not an entirely new 
concept. The rigid sectional templates often used for forming large curved surface 
structures also belong to the class of discrete surface tooling. Figure 1.11 shows one 
such tooling for forming stiffened panels used in an aerospace application. However, 
such a discrete tooling is not adjustable and can not be used for shapes other than the 
one for which it is built. 

The concept of adjustable discrete surface tooling was first given by Nakajima (Naka- 
jima, 1969). The concept was used to fabricate complicated dies for press work and 
electrodes for electrochemical machining. It consisted of number of thin wires bound 
together by a retainer into a bundle and the surface composed of each wire end was 
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used to simulate and replace the desired die surface (Fig. 1.12). Each of the wire is a 
surface element of an expected contour which can be pushed in the longitudinal direc- 
tion to a desired position by an ultrasonic vibrator. The positioning of each wire was 
carried out manually using a master model of the curved surface or by computer aided 
numerical control. Experiments were also carried out with a rubber sheet inserted 
between the blank and the die to avoid the scratching of the blank by wire elements. 
Constant diameter circular wires were used and the locking of wires was done with an 
external force applied perpendicular to longitudinal direction of the wires. The tooling 
is meant for small lot production or trial production and can be extended to die casting, 
compression molding and vacuum forming of plastics with appropriate modifications. 

Nakajima's work successfully demonstrates the use of discrete surface tooling for man- 
ufacture of sheet metal and plastic components. However, it does not deal with the 
detail study of tooling parameters in relation to the accuracy of the product shapes 
obtained. Once the pins forming the discrete die are positioned using a master model 
or numerical control, it is quite unlikely that the product with a desired shape can 
be obtained in the first trial of experiment. This is because the springback in terms 
of elastic recovery of the sheet after unloading manifests in all sheet metal forming 
processes. Moreover, the deformation of rubber sheet placed between the discrete die 
and the product being formed is likely to change the effective die shape at the time of 
forming. It requires a number of attempts by trial and error method in order to obtain 
products with a desired surface shape. 

An extensive research work on the concept of discrete surface tooling ha.s been carried 
out by Hardt et. al at the Laboratory for Manufacturing and Productivity at Mas- 
sachusetts institute of technology (Hardt and Gossard, 1980; Hardt and Webb, 1982; 
Webb, 1987; Webb and Hardt, 1991; Ousterhout, 1991; Hardt, 1992; Hardt et. al, 
1993). They propose the use of a variable configuration discrete element die, which is 
similar to Nakajima’s flexible die. Unlike the forming system proposed by Nakajima, 
they propose the use of discrete element die in a closed loop shape control system. The 
system forms a workpiece by an iterative procedure involving forming, measuring and 
die shape change until the shape convergence occurs. Figure 1.13 shows the conceptual 



Figure 1.12. Nakajima’s Discrete Surface Tooling. 
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sketch of their flexible sheet forming system. 

A rapidly reconfigurable discretized forming machine tool was developed on a large scale 
for this purpose. The machine tool consists of arrays of square pins with spherical ends, 
clamped in a housing by a hydraulic actuator with a force of 270 kN. The die shape 
can be programmed by a series of position servos that adjust a single die element at a 
time. The deformed shape of the part is measured after forming. The measured shape 
is compared with the desired shape to generate the shape error. The error is processed 
by a shape control algorithm to command a new geometry. The new die is then used 
to deform the workpiece to reduce the shape error. This forming cycle continues until 
an acceptably low shape error results. 

The reconfigurable discrete element die proposed by Nakajima and improved by Hardt 
el. al is primarily meant for forming sheet metal components. The work carried out 
by both these groups is based on the forming experiments using such a die, either in 
an open-loop or a closed-loop control system. Another approach to study the realiz- 
ability of shapes by a discrete element die is using numerical methods such as FEM. 
A computational study will be helpful in reducing the number of forming experiments 
required to arrive at the correct die shape. However, in the case of sheet metal forming 
processes a numerical study accounting for large strains, material properties of the 
sheet, springback and the tribological aspects associated with the process is difficult. 
In fact, the numerical simulation of sheet metal forming processes with continuous 
forming surfaces provided by the conventional dies is still an active area of research 
(Cho, 1993; Hasan, 1993; Zhou, 1993). 

In the case of composite fabrication processes such as layup, the complexities associated 
with the sheet metal forming are absent. In this case the forming forces are also 
considerably less. However, the thermal effects are likely to be present in all composite 
fabrication processes. But, in the case of wet or prepreg layup with room temperature 
curing these effects can be neglected. This means that, once the tool surface in terms of 
deformed flexible sheet is realized, the realization of the part is ensured. It implies that 
the product shape here is completely defined by the mold shape and the realizability 
of the shape, is independent of the process. These simplifications make the numerical 
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Open Loop Sheet Forming System (Nakajima) 



Closed Loop Sheet Forming System (Hardt et al.) 


Figure 1.13. Open and Closed Loop Shape Control Systems for Sheet Metal 
Forming with Variable Configuration Die. 
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study of the flexible surface tooling feasible for applications like composite layup. 


1.8 Objectives and Scope of Work 

The overall objective of the present work is to carry out the numerical study of the 
flexible surface tooling concept in order to evaluate its applicability to composite layup 
process in a prototyping environment. The specific objective of this research can be 
stated as follows: 

• To identify the shape conforming processes. 

A class of manufacturing processes have been identified where the shape realization is 
due to geometric transformation from an original shape to a tool shape. 

• To propose the discrete surface tooling concept to various shape conforming pro- 
cesses. 

A discrete surface tooling which is based on the concept of discrete approximation to 
a continuous surface, has been extended for shape conforming processes in general and 
composite manufacturing processes in particular. 

• To study the realizability of shapes by flexible surface tooling. 

A numerical study of the realizability of simple and complex shapes by the proposed 
flexible surface tool has been carried out. A computer program which deals with 
the deformation analysis of nonlinear elastic membranes in multiple contact has been 
written and used for this purpose. 

• To study the effect of tooling parameters on the surface shape accuracies. 

The effect geometric and material properties of the flexible sheet on the shape accura- 
cies has been studied with an objective to identify the range of these parameter for a 
given application. 
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• To identify and study the parameters associated with discrete approximation to 
a continuous surface. 

The tooling parameters related to surface discretization such as the number and the 
layout of indentors, the section and the contact geometries of the indentors play a 
major role in obtaining accurate surface shapes. The effect of these parameters on the 
accuracy of surface shapes has been studied with an objective to identify the suitable 
range of theses values for fabricating a flexible surface tooling machine. 

• To study the static shape optimization of the flexible sheet in the proposed tooling. 

The static shape of the deformed flexible sheet in the proposed tooling is used as 
a base surface for tooling applications. Apart from the other parameters associated 
with the tooling, the deformed shape of the flexible sheet depends on the heights by 
which indentors are raised. The sensitivity of the static shape to the indentor heights 
is an important aspect. This has been dealt with an objective to improve the shape 
accuracies obtainable. 

• To study the realizability of shapes having negative Gaussian curvature. 

The proposed scheme of tooling ha^ been modified for realizing the shapes having 
negative Gaussian curvature. A numerical study of this scheme of tooling has been 
carried out in order to evaluate the realizability of shapes. 

The concept of discrete surface tooling has been extended to some of the shape con- 
forming processes as a part of the present work. However, the present numerical study 
of the proposed tooling is restricted to only those processes where a single finished sur- 
face can be obtained. Another assumption of the present work is that the realizability 
of a component shape is assumed to be dependent on the the realization of tool surface 
and is assumed to be independent of the process. Such an assumption is valid only 
when the forming forces are considerably less and the thermal effects are insignificant. 
This is true with processes like composite layup with room temperature curing. 
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1.9 Organization of the Work 

A computational study of the proposed flexible surface tooling, for realizing shapes 
having positive and non-positive Gaussian curvature is the main aim of the present 
research work. This will help in knowing the shapes realizable, identifying the tooling 
and discretization parameters for obtaining products with improved shape accuracies. 

In the proposed tooling system, deformed configuration of the flexible elastomeric sheet 
is used as the mold surface. Hence, the deformation analysis of the flexible sheet forms 
the central idea of the present work. This involves solving a finite elasto-static contact 
problem with multiple contacts. Chapter 2 presents the formulation and the numerical 
scheme used for solving this problem. The chapter describes the modeling of the 
flexible elastomeric sheet for the purpose of analysis, finite element formulation of the 
problem, the computational procedure and the contact searching scheme. The chapter 
also presents the validation of the computer code developed for the purpose of analysis. 

Chapter 3 starts with the description of the steps involved in the shape realization pro- 
cedure followed here. This consists of a preprocessor to determine the indentor heights, 
the deformation analysis and the post processing for the evaluation of error between 
the desired and realized surface. The chapter presents the results of the numerical 
experiments in terms of realization of singly and doubly curved surface shapes. The 
results of the study of various parameters associated with the flexible sheet such as 
geometric and material properties of the sheet, sheet dimensions etc. on the accuracy 
of shapes realized are presented. The effect of external loading on the realized shapes 
is also dealt as a part of this study. 

A discrete approximation to a continuous surface is an important ingredient in the 
proposed tooling. Chapter 4 is concerned with the study of parameters associated 
with the surface discretization. These include the number of indentors, the layodt 
of indentors, the section and contact geometries of the indentors etc. The objective 
of this study is twofold. The main purpose of this study is to know the effect of 
discretization parameters in order to improve the accuracy of realizable shapes. The 
secondary objective of such a study is to identify the range of discretization parameters 
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for a surface tooling machine to be fabricated. 

Chapter 5 of the thesis deals with the static shape optimization of the flexible sheet 
for the indentor heights. In the absence of the flexible elastomeric sheet between the 
mold surface and the part, the problem of discrete approximation to a given continuous 
surface turns out to be a trivial one. But the presence of flexible sheet on the discrete 
surface of the mold changes the effective mold shape. This change is significant for 
higher values of the thickness of the flexible sheet. In such cases, the sensitivity of the 
static shape of the flexible sheet for indentor heights forms an important and interesting 
aspect of the present study. This aspect has been investigated with an objective to 
minimize the error between the desired and realized surface. 

A flexible surface tool in which the flexible sheet is constrained by some or all the 
indentors bonded to it has been proposed as another version of the tooling. This 
version of the tooling can be used for realizing shapes having both synclastic and 
anticlastic regions. A numerical study of this version of tooling forms the contents of 
Chapter 6. The main objective of this study is to evaluate the realizability of shapes 
having negative Gaussian curvature. A secondary objective of this study is to compare 
this version of tooling with the earlier version dealt as a part of Chapter 3. 

The conclusions of this work are presented in Chapter 7. The future scope of the work 
in terms of analysis, design and fabrication of the proposed tooling are discussed in 
detail. The possibility of extending the concepts and results which are the outcome of 
the present work to applications in other areas of engineering is also discussed. 



Chapter 2 


DEFORMATION ANALYSIS OF 
FLEXIBLE SURFACE TOOL 


2.1 Introduction 

A flexible surface tooling system has been proposed in the previous chapter for shape 
conforming processes. In the proposed tooling system, the deformed configuration of 
the flexible elastomeric sheet is used as a base surface for the tooling requirements. The 
elastomeric sheet undergoes deformation due to indentation by a matrix of indentors. 
Knowing the deformed configuration of the flexible sheet for a given configuration of 
the indentors is an important aspect in the above mentioned tooling. If the thickness 
of the flexible sheet is small, and if the change in sheet thickness due to indentation 
is uniform, and also if the sheet is in contact with all the indentors then the deformed 
configuration of the flexible sheet can be known by geometric methods. But, this is 
not true in most of the cases. The complexities associated with the deformation of 
rubber-like materials make this problem a non trivial-one. 

A mechanics-based approach accounting for the large deformations of rubber-like mate- 
rials and the contact boundary conditions is essential to solve the above problem. The 
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presence of contaci constraints together with nonlinear elastic nature of the flexible 
sheet makes the solution of this problem possible only by numerical methods such as 
FEM. This chapter describes the formulation of the above elasto-static contact problem 
and the numerical scheme used for solving it. 


2.2 Scope of the Analysis 

The flexible sheet used in the proposed tooling is subjected to a considerable stretching 
due to indentation. The sheet is made of rubber-like material and can be modeled as 
an isotropic hyperelastic medium for the purpose of analysis (Oden, 1972). Further, 
most of the rubber-like materials are treated as incompressible because of their capa- 
bility of sustaining severe deformations without an appreciable change in volume. The 
flexible sheet here is of uniform thickness before deformation and it’s thickness is small 
compared to it’s other dimensions (less than 1%). As the state of stress in the sheet is 
essentially two dimensional in this case, the sheet can be idealized as a thin membrane 
with no bending effects. Thus, in the present analysis the flexible sheet is treated as a 
thin incompressible isotropic hyperelastic membrane. 

The present computational analysis of the flexible sheet is limited to calculation of dis- 
placements only, and no attempt has been made to evaluate the stresses. It is assumed 
that the flexible sheet is capable of taking the applied strains without any failure. The 
contact friction between the flexible sheet and the indentors is not considered as a part 
of the present analysis. The flexible sheet is subjected to indentation on one side and 
the other side is acted upon by some external forces due to material processing. These 
external forces are small in magnitude for processes like layup. Secondly, the surface 
finish of the indentors is high. These factors contribute in reducing the contact friction 
between the sheet and the indentors. 
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2.3 Material Description 

The deformation process involving hyperelastic materials is reversible and isothermal 
in nature. These materials are characterized by the presence of an elastic potential, w 
known as strain energy density which represents the strain energy per unit volume of 
the undeformed body. For rubber-like materials which can be considered as isotropic, 
the strain energy density is of the following form (Oden, 1972). 

w = w{Iul 2 ,h) (2.1) 

where /j, I 2 and I 3 are the three strain invariants. The incompressibility condition can 
be imposed by treating the third strain invariant I 3 as unity. Accordingly the strain 
energy function for incompressible material reduces to following form. 

w = w{Ii,l 2 )-, h = l. (2.2) 

Various approximations of w for specific materials have been proposed. A general form 
of these expressions is given by, 

00 00 

» = E E - 3)'(h - 3)'; h = 1. (2.3) 

t=0 j=0 

The most widely used special case of (1.3) is the Mooney-Rivlin material description 
(Mooney, 1940; Rivlin, 1948). It is used for characterizing rubber-like materials under- 
going large strains. It is given by, 

u, = C, (/a - 3) -h (72 (/2 - 3) ; h = l. (2.4) 

where (7i, and C 2 are the material constants. For C 2 = 0, the strain energy function 
w, reduces to a standard neo-Hookean form given by. 


in = C(/i-3) 


(2.5) 
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2.4 Earlier Work 

The first significant theoretical study of the large deformations of rubber-like materials 
was carried out by Adkins and Rivlin (1952) using the neo-Hookean and Mooney forms 
of strain energy function. More details of this work can be found in the book by Greens 
and Adkins (1960). Subsequent contributions in this direction in terms of analytical 
studies were made by Klingbeil and Shield (1964), Hart-Smith (1966), Foster (1967) 
and Hart-smith and Crisp (1967) etc. Finite element methods for the analysis of finite 
deformations of elastic bodies were first developed by Becker (1966) and Oden (Oden, 
1966; Oden, 1967; Oden and Sato, 1967) independently. They attempted the plane 
stress problem involving large deformations of rubber-like materials. Subsequently, the 
problems of finite plane strain, axisymmetric and three dimensional cases were solved 
by finite element method (Oden, 1968; Oden and Key, 1970; Oden and Key, 1971). 
Much of this early work in this direction can be found in the books of Oden (1972, 
1984). 

Contact problems involving elastic bodies have attracted considerable attention since 
the work of Hertz (1882) in the last century. Problems of solid and structural mechanics 
with contact constraints have in recent years been the subject of a number of both 
theoretical and numerical investigations. There is enough literature, now available on 
contact problems and the more details of which can be found in the texts/monographs 
by Gladwell (1980), Kikuchi k Oden (1988) and Zhong (1993). As the present work 
deals with a contact problem involving hyperelastic membrane, only those references 
which are of direct 'relevance to the present problem are discussed in more detail. 

One of the first contact problems involving elastic membranes was attempted by Yang 
and Hsu (1971). They solved the problem of axisymmetric indentation of a circular 
membrane by a rigid sphere. The problem was solved by dividing the deformed mem- 
brane into two regions: a region of contact with the sphere and a region which is free 
of external load except boundaries. The common boundary between the two regions 
is not known a priori. The nonlinear membrane equations are applied specifically to 
each of the regions satisfying the continuity of stress and deformation at the common 
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boundary. 

The inflation analysis of highly elastic membranes in the presence of rigid constraints 
has been carried out by Feng et al. (Feng and Yang, 1973; Feng, 1974; Feng, 1975). 
The coordinate functions that describe the deformed configurations of the membrane 
are assumed to be represented by a series of geometrically admissible functions with 
unknown coefficients. The shape and the boundary of the contact region and the 
configuration of the deformed membrane under both inflation and indentation were 
found by minimizing the potential energy. The optimization process is carried out by 
using the Fletcher-Powell iterative descent method subjected to inequality constraint 
conditions. 

The confined axisymmetric inflation of a hyperelastic cylindrical membrane with fixed 
ends is examined by Khayat et al. (1993). Assuming a Mooney-Rivlin constitutive 
model the problem is formulated to account single and multiple contacts. The tan- 
gency conditions between elastic membrane and solid boundary is respected but the 
formulation does not account for large deflections. 

In the context of thermoforming process, the finite element analysis of rubber-like 
membranes with contact constraint has been dealt by a few investigators (Charrier et 
al., 1987; deLorenzi and Nied, 1987; Charrier et al., 1989). Thermoforming process 
consist of inflating or pressing a heated polymer membrane against a relatively cold 
mold. In these works the hot polymer is modeled as a nonlinear elastic incompressible 
membrane. The contact between polymer and the mold is treated as frictionless and 
non-slipping. The work of Charrier and associates deals with the deformation analysis 
of axisymmetric and non-axisymmetric shapes using one dimensional and three noded 
triangular finite elements respectively (Charrier et ah, 1987; Charrier et ah, 1989). 
De Lorenzi and Nied (1987) use finite element program ADINA for modeling thermo- 
forming and blow molding. More details of the application of finite elements in the 
thermofonning process can be found in Zamani et al. (1989). 

Ogden’s strain energy function based on principal stretches has been used by Gruttman 
and Taylor (1992) for finite element formulation of contact problem involving elastic 
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membranes. Inflation analysis of a torus constrained between two plates has been 
carried out by them using Newton’s algorithm and penalty based contact constraint. 

In the contact problems discussed above the rubber-like bodies are in the membrane 
or plane state of stress. The contact problems of hyperelastic materials in the state of 
finite plane strain, finite axisymmetric deformations and the three dimensional state of 
stress have been attempted by many investigators. The notable among them are the 
works of Kikuchi & Song (1980), Endo et al. (1984), Batra (1985), Sussman & Bathe 
(1987), Bass (1987), Haggblad k Nordgren (1989) etc. 


2.5 The Present Methodology 

Most of the problems dealing with rubber-like materials in contact have been attempted 
by finite element method. The incremental-iterative strategies based on total or up- 
dated Lagrangian concepts are commonly used to solve this nonlinear problem. The 
material incompressibility and contact constraints are dealt either by mixed formula- 
tion or penalty based approach (Gadala, 1986). However, when the problem involves 
frictionless contacts and plane or membrane state of stress it can also be attempted by 
a mathematical programming approach. This approach, is commonly known as energy 
search, and consists of searching for the equilibrium configuration of the body which 
corresponds to a state when the potential energy is minimum. Such an approach is 
widely used for solving linear elasto-static contact problems (Chand et al., 1976; Hung 
and Saxce, 1980; Kamat and Hayduk, 1980; Westbrook, 1982; Sayegh and Tso, 1986; 
Johnson and Quigley, 1989; Lee and Kwak, 1989). The mathematical programming 
approach has been followed by a number of investigators for the deformation analysis 
of rubber-like membranes with or without contact constraints (Levinson, 1965; Becker, 
1966; Feng and Yang, 1973; Feng et al., 1974; Tielking and Feng, 1974; Feng and 
Huang, 1974; Feng and Huang, 1975; Rigbi and Hiram, 1981). While, Becker uses the 
finite element discretization together with the energy search, other investigators follow 
Rayleigh- Ritz type of procedure with continuous trial functions. 

In the present work the problem of deformation analysis of flexible sheet has been 
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attempted by finite element discretization and potential energy minimization in the 
presence of contact constraints. The approach is same as that of Becker (1966) except 
for the contact constraints. Three noded triangular elements are used for the finite 
element discretization. Mooney-Rivlin model is used for the strain energy calculations 
with incompressibility condition. The incompressibility condition has been dealt by 
analytical treatment which is possible in the case of plane or membrane state of stress 
(Gadala, 1986; Oden and Sato, 1967). The conjugate gradient method with Polak- 
Ribiere correction is used for energy minimization. The contact constraints are treated 
using a penalty approach. 


2.6 Theoretical Development 

In this section, kinematics of thin membrane geometry, membrane strain tensor and 
potential energy expressions are presented. Lagrangian description of motion referred 
to a set of convected coordinates is used to describe the membrane kinematics. The 
theoretical formulation used here is similar to that of Oden & Sato (1967). As the 
analysis is restricted to calculation of displacements only, an explicit computation of 
hydrostatic pressure is not necessary. Because, in an incompressible deformation pro- 
cess the hydrostatic pressure component affects only the stresses but not strains or 
displacements. 

Consider the geometry of deformation of a thin homogeneous elastic membrane of 
uniform thickness ho. The un deformed and deformed configurations are represented 
by Oo and f2 respectively (Fig. 2.1). The region fio is bounded by the planes iY3±ho/2, 
the middle surface being X 3 = 0. The position vectors x = xjii -t- X 2 i 2 + 2:313 and 
X = JAiii -t- X 2 i 2 4- -Asia of the deformed and initial configurations respectively, are 
related through expression given by. 


X = X + u 


( 2 . 6 ) 


where u is the displacement vector. 
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The Lagrangian strain tensor is now given by the relation 

1 / dxk dxk 


— o 


2\dXidXi 


( 2 . 7 ) 


where is the Kronecker delta. Thickness of the membrane after deformation, h, is 
given by 


h = \ho ( 2 . 8 ) 

where A is the extension ratio at the middle surface in a direction normal to the middle 
surface. Assuming the membrane to be thin, strain tensor can be written as 

if dua , dup ^ dui dui ^ 

~ 2 ^ ^ j 

7i3 = 0 723 = 0 733 = ^ (A^ - l) (2.9) 

where a, /? = 1,2 and i = 1,2,3. The three strain invariants /i, 72 , h are given by 

7i = A^ + 2(l+7aa) 

h = —A'* + A^/i + 

h = A^ (1 + 27a„ + 2ea/5eA;u7QA7/?M) ( 2 - 10 ) 


where a, ^5, 7 , /z = 1,2 and Cap and e^ti are the two dimensional permutational symbols 
(ei 2 = 1, €21 = —1, eii = 622 = 0). Imposing the incompressibility constraint 73 = 1 , 
the above expressions can be written as, 

h = A^ + 2(l+7aa) 

I 2 = — A^ + A^7i + 

X‘‘ 

A — (I "h 270101 T 2eot^eAfj7ofA7)3A‘) 


( 2 . 11 ) 
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The potential energy for the deformed membrane, 11 is given by 

U = U -W = [ w{hJ 2 )dv- [ t-\ds (2.12) 

Jv J S 

where U is the strain energy of the deformed membrane and W the virtual work of 
external traction forces, t. v is the kinematically admissible displacement field and s 
is the part of boundary on which tractions are specified. 


2.7 Finite Element Formulation 

The finite element idealization used here consists of replacing the middle surface of the 
membrane by a patchwork of plane triangular elements and loading it with equivalent 
concentrated loads at the nodes. Local state of stress here is considered to be that 
of plane stress (as = 0). The analysis further assumes that the principal stresses are 
tensile everywhere {cri >0, <72 > 0). This is the central idea on which the works of 
Oden and Sato (1967) and Charrier et al. (1989) are based. 

Consider a triangular element with nodes A, B and C in a coordinate space, defined 
by a triad of unit vectors (ii 12 is). The element undergoes deformation to take a new 
shape U® with nodes a, h and c (Fig. 2.2). In order to calculate the strain energy 
due to stretching of the membrane, it is convenient to compare the two configurations 
in a common plane. For this purpose, undeformed element is placed in the plane of 
deformed one. Such a transformation does not affect the strain energy calculations as 
it involves only rigid body translations and rotations. 

We define a local coordinate system for the undeformed element by a triad of unit 
vectors {E}^ = (El E2 Ea). The unit vector Ei is directed from node A to B. The 
unit vector E3 is perpendicular to the plane of the element and E2 is the cross product 
of E3 and El- Similarly, we attach a local coordinate system for the deformed element 
which is defined by the triad of vectors {e}^ = (ei 62 63). The two local coordinate 
systems are related to global coordinate system {i} = (ii 12 13) by transformation 

matrices [R] and [r] given by, 




Figure 2.2. Deformation of a Triangular Membrane Element in Space. 
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|«l = {E}{i)^ M = {a}{i}’' (2.13) 

Figure 2.3 shows the undeformed element (after transformation) and the deformed 
element in the local coordinate system of the latter. Let {u}^ = {ua ub wc} and 
{u}^ = {va vb vc} are respectively the nodal displacements of the triangular element 
in ei and 62 directions respectively. The displacements u and u of a point within the 
element, using a linear variation is given by, 

u = Naua + Nbub + Ncuc = {u}'^ [N] 

V = Nava + Nbvb + Ncvc = [N] (2-14) 


where Na, Nb and Nc are shape functions given by. 


NA{x,y) 

Nc{x,y) 


1 

2A 

1 

2A 

1 

2A 


[{yc - vb) (xb -x)- {xc - xb) [yB - y)] 
{[Va - yc) {xc -x)- {XA - Xc) {yc - y)] 
[{yB - yA) {xA -x)- {xb - xa) {yA - y)] 


2A = {xc - xb) {yA - yc) - {xa - xb) {yc - yB) 


(2.15) 


Using Equation (1.15), one can derive the the following. 


1 “ = -;h'\.{yB~yc)uA + {yc-yA)uB + {yA-yB)uc] 

ax zZa 

fx ~ ^ [(ys - yc) Va + {yc - yA) vb + {yA - yB) vc] 

^ = :^[{xc-xb)ua-+{xa-xc)ub + {xb-xa)uc] 

= ■^[{xc-xb)va + {xa-xc)vb + {xb-xa)vc] 
ay zZX 


( 2 . 16 ) 


the strain energy, U can now be evaluated by substituting the results of (1.16) in the 



Figure 2.3. Comparison of Deformed and Undeformed Elements in the 
Plane of Deformed Element. 
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expressions of membrane strain tensor given by Equation (1.9) and using Equations 
(1.11) and (1.12). 

The external forces in this problem are in the form of pressure loads acting on the 
surface element. As the membrane deforms not only does the direction of this pressure 
change but also the area on which it acts changes. To account this, it is first assumed 
that the dimensions of each element are sufficiently small that the pressure is essentially 
unform over the surface of individual elements. It can be shown that consistent with 
the linear variation of displacement assumed within an element and the positive sense 
of pressure, the equivalent out of plane forces at the element nodes are 

=F,b^ (2.17) 

where p is the applied pressure and is the deformed element area. The global 
components of the nodal forces can now be obtained by the following transformation. 



( 2 . 18 ) 


2.8 Contact Constraints 

The deformed configuration of the flexible sheet due to indentation is such that the 
geometric constraints are not violated. In other words any solution to the equilibrium 
configuration of the membrane is such that no interpenetration of the indentor and 
membrane materials is allowed. This condition on the kinematically admissible dis- 
placement field of the flexible sheet can be imposed by penalty approach. The penalty 
being proportional to the amount of violation. 

The above contact conditions which are not known a priori needs to be evaluated by 
contact searching. The two models which are often used for evaluation of contacts 
are node-to-node interface model and node-to-segment interface model (Zhong, 1993). 
The node-to-segment contact searching scheme is generally used for bodies undergoing 
large displacements and rotations which consist of discretizing the contacting body in 
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terms of finite elements and also discretizing the target body in terms of polygons. 
Then, each of the potential contacting nodes of the contacting body are checked with 
each of the segments of the target body. 

A slightly modified version of the node-to-segment algorithm has been used here which 
consists of, checking for each of the potential contacting nodes of the contacting body 
with the surface of the target body without any discretization. The target bodies here, 
namely, the indentors are rigid and axisymmetric. They can be defined in terms of 
surface equations and this definition of surface is directly used for contact evaluation. 
A number of minimax tests are used to identify potential contacting nodes and to avoid 
any unnecessary computation involved. 


2.9 Computational Procedure 

The computational procedure used here consists of, minimizing the total potential 
energy of the system subjected to given geometric constraints. For this purpose, the 
domain Oq is discretized into three noded triangular finite elements. The total number 
of nodes and elements being JV and E respectively. The number of variables in the 
minimization problem are 3N'; N' being the total number of unconstrained nodes. 
The objective function is the total potential energy of the system summed over E 
elements. Potential energy of the individual element consists of the strain energy due 
to deformation obtained from undeformed and deformed configurations of the element 
after transformation and the potential of loads, if any, on the element. 

To solve this constrained optimization problem efficiently, the sequential unconstrained 
minimization technique (SUMT) has been used. SUMT transforms the constrained 
optimization problem to an unconstrained one by penalizing the constraint violation 
if any. The conjugate gradient algorithm with Polak-Ribiere correction is adopted as 
the minimization technique. This technique has been used by other investigators also 
for a successful solution of contact problems. (May, 1986; Dilintas et ah, 1988). The 
method of multipliers (MOM) type of penalty which is equally effective in both feasible 
and non-feasible solution domains is used for transforming constrained optimization 
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problem into an unconstrained one (Schuldt et al., 1977; Reklietis, 1989). The computer 
program developed for this purpose is written in C and is implemented on DEC-3000 
platform. 

It is convenient to assume the initial or undeformed configuration as a starting solution 
to the above optimization problem. A starting solution close to the equilibrium con- 
figuration is found to save a lot of computational time. The analysis has been carried 
out with different starting solutions and in every case the solution converged to the 
equilibrium configuration. These experiments justify the uniqueness of the solution in 
this case. The experiments were also carried out with single and multiple load steps. 
Based on the results of these experiments, it has been found that the convergence of 
the solution in this case is independent of the load step. This further strengthens the 
uniqueness of the solution and the smooth convergence of the optimization process in 
this case. 


2.10 Validation 

The correctness of the formulation and the numerical scheme used here has been veri- 
fied by running the computer program developed, for three standard test cases taken 
from the works of Becker (1966), Feng k, Huang (1974) and Feng et al. (1974). The 
first test case deals with the problem of in-plane stretching of an elastic sheet. A thin 
square membrane is fixed along one of its edges and the opposite edge is moved to 
stretch the membrane to twice its original length. The purpose of this example is to 
verify the correctness of strain energy calculations and the convergence of the solu- 
tion process. The ratio of first and second Mooney constants, C\IC 2 is taken as 8.0 
for the analysis. Figure 2.4 shows the deformed profile of the stretched membrane. 
The deformed geometry has been compared with the results of Becker (1966) and the 
maximum deviation in displacement has been found to be the order of 0.8% (Fig. 2.5). 

The second test case is concerned with the inflation of a square membrane for various 
applied pressures. The results are compared with the work of Feng & Huang (1974). 
Figure 2.6 shows the deformed configuration of the membrane at various inflation 
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Figure 2.4. In-plane Stretching of a Thin Elastic Membrane. 





Figure 2.5. In-plane Stretching of a Thin Elastic Membrane: A comparison 
with Becker’s Solution. 
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pressures. The inflation pressure values are measured in terms of non-dimensional 
pressure value P* given by, 



where p, Ci and h are the applied pressure, the non-dimensional Mooney constant and 
the thickness of the membrane respectively as defined in the Feng & Huang (1974). 
The square membrane is of the size 2a. The ratio of first and second Mooney constants, 
Ci/Cj is taken as 5.0. Figure 2.7 shows the variation of the non-dimensional height 
of the middle point of the membrane, Fl for different non-dimensional pressure value, 
P*. The results of the work carried out by Feng & Huang are also shown in the figure 
by means of small circles. 

The first two test cases conform the correctness of strain energy calculations and the 
potential of loads. The third test case is an extension of the second one where the 
inflation of a square membrane is carried out in the presence of a rigid geometric 
constraint. The ratio of first and second Mooney constants, CijC 2 is taken as 10.0 
here. The geometric constraint is in the form of a rigid flat surface placed above the 
membrane. As the applied inflation pressure increases, the membrane comes in contact 
with the flat constraint which acts as a barrier. The deformed profile of the center-line 
of the membrane is plotted in Figure 2.8. The results are compared with that of Feng 
and associates (1974). The deformed profile of Feng et al. are indicated by small 
circles. 



Figure 2.6. Inflation of a Square Membrane. 




Y3 Non Dimensional Center Point Height 
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Figure 2.7. Inflation of a Square Membrane for Various Applied Pressures. 
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for the indentors is assumed which is also a datum for measuring the heights of the 
individual indentors. In their base position the indentors do not make any contact with 
the flexible sheet. As the indentors are raised above their base position they come in 
contact with the flexible sheet and deform it. 

The main objective in the proposed tooling system is to deform the flexible sheet to 
a desired shape. This involves, knowing the heights by which the individual indentors 
are to be raised in order to deform the flexible sheet to a required shape. There can be 
two approaches to attempt this problem. In the first approach, the heights by which 
the indentors are to be raised to realize a given surface are known in advance. This is 
an open loop approach. Using a matrix of indentors a discrete approximation to a given 
continuous surface can be built by simple geometric methods. The relative heights of 
the indentors which form such an approximation can be used to realize the desired 
surface in the proposed tooling. 

In the second approach, the heights by which the indentors are to be raised are not 
known in advance. The process of shape realization starts with some initial assumed 
heights of the indentors. The deformed shape of the flexible sheet is measured for this 
set of heights and the error between the desired and realized surfaces is fed to some 
shape control algorithm which suggests new positions of the indentors. This procedure 
is repeated until the shape error ceases to improve beyond certain limit. Such an 
approach is known as closed loop approach. The numerical experiments carried out as 
a part of this chapter are based on the first approach. Some of the results obtained by 
following the closed loop approach are presented in the subsequent chapters. 


3.2 Shape Realization Software 

In order to study the shape realization by the proposed tooling, a set of computer 

programs have been written. The software developed for this purpose is in three 

\ 

modules: a preprocessor, a deformation analysis program and a post-processor. The 
input to the shape realization software consist of the following. 
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• ,The shape to be realized or the desired surface shape 

• Dimensions of the flexible sheet 

• Sheet thickness 

• Material properties of the sheet 

• Geometry of the individual indentors 

• Layout or the arrangement of indentors 

Knowing the heights by which indentors are to be raised to approximate a given contin- 
uous surface, is the first step in the process of shape realization. The preprocessor part 
of the software deals with this aspect. For this purpose, a common coordinate system 
is used to define the desired shape and the indentor positions. Each of the indentors 

is now raised incrementally till it makes a contact with the desired surface. After all 

/ 

the indentors are raised in a similar manner, it forms a discrete approximation to the 
given continuous surface. 

The heights of the indentors decided as part of the preprocessor together with various 
tooling parameters mentioned above, form the input to analysis module. The deforma- 
tion analysis program is based on the theoretical formulation and the computational 
procedure described in the previous chapter. The analysis program uses the middle 
surface of the flexible sheet as a reference surface, whereas, the upper surface of the 
flexible sheet is the desired one. This offset can be easily accounted using Equation 
(2.8) as the deformation of the flexible sheet proceeds. 

As a part of the post-processor, the deformed configuration of the flexible sheet obtained 
as a result of analysis is compared with the desired surface shape to evaluate the shape 
error. The realizability of a given shape here is measured in terms of the root mean 
square (RMS) error and the average absolute (AVAR) error between the desired and 
realized surface shapes. Before such a comparison is made, one of the surface shapes 
needs to be localized with respect to the other. Localization refers to the problem of 
determining the optimum position of one surface with respect to other for a minimum 
distance norm. The localization procedure followed here consists of positioning the 
deformed surface shape with respect to desired one by exhaustive search. 
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3.3 Numerical Examples 

Five examples of shape realization are presented in this section. Three of which are 
concerned with the surface patches which are doubly curved. The first example deals 
with the realization of a simple cylindrical patch. The parametric representation of the 
patch is given by, 


X = 2.8625 cos{u) 
y = 1.0*v — 0.5 

z = 2.8625 sm(u) - 2.7625 (3.1) 

where 1.395 < u < 1.746 and 0.0 < u < 1.0 

In order to realize the above surface, a rectangular sheet of size 1.6 m by 1.2 m and 
thickness 4 mm is indented with a total of 121 indentors. Two of the edges of the 
flexible sheet of 1.2 m length are constrained. The square indentors of 99 mm side 
with spherical end (radius=100 mm) are used. The indentors are placed in an area of 
1.1 m by 1.1 m forming 11 rows, each row consisting of 11 equally spaced indentors 
(Fig. -3.1). The ratio of first and second Mooney constants, C 1 /C 2 is taken as 4.0. The 
analysis has been carried out with a total of 3072 triangular elements (Fig. 3.2). 

Figure 3.3 shows the desired and deformed surface patches. The figure shows only 
that region of the deformed surface which is in contact with the indentors. Figure 
also depicts the error between the desired and deformed surface patches at higher 
magnification. 

In order to realize the above surface patch, the analysis is also carried using a square 
sheet of size 1.6 m with all the four sides constrained. The average absolute (AVAB) 
error and the root mean square (RMS) error between the desired and deformed surfaces 
for the two cases are indicated in Table 3.1. 
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Figure 3.2. Finite Element Discretization of the Membrane (Numerical 
Example I). 
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Boundary 

Average Absolute 

Root Mean Square 

Condition 

Error 

Error 

Two Sides Constrained 

0.2595 mm 

0.3248 mm 

Four Sides Constrained 

0.3261 mm 

0.3931 mm 


Table 3.1 AVAB and RMS Errors (Numerical Example I). 

In the second example a rectangular membrane of thickness 5.0 mm is indented with 
multiple indentors to simulate an elliptic cylindrical surface patch. Circular indentors 
of radius 49.5 mm are used (Fig 3.4). The other tooling parameters being same as 
that of previous example. The parametric equation of the desired surface is given as 
follows. 


X — 1.51625 cos(u) 
y = 1.0*u — 0.5 

2 = 0.758125 sm(u) - 0.658125 (3.2) 

where 1.0534 <u< 2.0882 and 0.0 < u < 1.0 

The desired and deformed configurations of the membrane are shown in Figure 3.5. 
The AVAB error and the RMS error between the desired and deformed surfaces are 
given in Table 3.2. 


Boundary 

Average Absolute 

Root Mean Square 

Condition 

Error 

Error 

Two Sides Constrained 

0.6936 mm 

0.8611 mm 

Four Sides Constrained 

0.7591 mm 

0.9350 mm 


Table 3.2 AVAB and RMS Errors (Numerical Example II). 


The third numerical example deals with the realization of a doubly curved surface 
patch. The equation of the surface patch expressed in an explicit form is given as, 


z = 0.08 [(41/3)* - l]'‘ [(4!,/3)* - l]"* 


(3.3) 
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where -0.5 < a; < 0.5 and -0.5 < j/ < 0.5 

In order to realize this surface patch, a square membrane of size 1.6 m by 1.6 m and 
thickness 5 mm is constrained at all the four edges and indented with a total of 121 
indentors. Other input data is same as that of example I. The desired surface, deformed 
configuration and the error profile are shown in Figure 3.6. The AVAB and RMS errors 
are given in Table 3.3. 


Average Absolute 

Root Mean Square 

Error 

Error 

0.3187 mm 

0.3975 mm 


Table 3.3 AVAB and RMS Errors (Numerical Example III). 


In the next numerical example, an asymmetrical bicubic Bezier patch has been realized. 
The characteristic polyhedron which defines the Bezier patch is given in Table 3.4. 
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Table 3.4 Bezier Characteristic Polyhedron (Numerical Example IV). 
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A square ineml)ranr with side 1.8 m has been used for the realization of this patch and 
the indentors are placed in an area of 1.1 m by 1.1 m. The results are presented in 
Figure 3.7 and the AVAB and RMS errors are given in Table 3.5. 


Average Absolute 
Error 

Root Mean Square 
Error 

0.2041 inm 

0.2531 mm 


Table 3.5 AVAB and RMS Errors (Numerical Example TV). 

Lastly, the results of realization of another doubly curved surface patch are presented 
in Figure 3.8. The curvatures of the surface patch here are much higher than those of 
the previous examples. The equation of the surface patch expressed in an explicit form 
is given as, 


2 = 0.25 [l-(5x/4)^] [l-(5y/4)''] (3.4) 

w'here —0.5 < x < 0.5 and —0.5 < y < 0.5 

The AVAB and RMS errors between the desired and the realized surface patches after 
localization are given in Table 3.6. 


Average Absolute 

Root Mean Square 

Error 

Error 

1.1949 mm 

1.4755 mm 


Table 3.6 AVAB and RMS Errors (Numerical Example V). 

The values of the average absolute and root mean square errors presented above are only 
typical values and not necessarily the best values obtainable. The error between the 
desired and deformed surface patches can be further reduced by selecting the optimum 






Figure 3.8. Realization of a Doubly Curved Surface Patch (Numerical Ex- 
ample V). 
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valiu'K of the parameters. The subsequent sections and chapters deal with this 

aspect in detail. 


3.4 Effect of External Loading on the Shape Ac- 
curacy 

The numerical examples presented in the previous section assume that there is no 
external loading on the flexible sheet. The deformation of the flexible sheet in these 
cases is only due to indentation. In practice, some form of an external loading on the 
flexible sheet can be expected. This may be due to loads encountered as a result of 
the material processing. The shape realization in the presence of external loading then 
becomes an important aspect. This aspect has been dealt as a part of the present 
section. 

The numerical experiments were carried out for two forms of loading acting on the 
flexible sheet. These are uniformly distributed and dead loading. With the uniformly 
distributed loading (UDL), the direction of the loading changes with the deformation 
of the. flexible sheet. The loads are assumed to be acting on the flexible sheet in an 
area where the indentors are placed. The effect of these two types of loads has been 
studied for a variety of surface shapes. Figures 3.9 to 3.12 present the results of the 
same for a typical surface shape. The thickness of the flexible sheet used here is 4 mm. 

Figure 3.9 shows the variation of average absolute (AVAB) error with an increasing 
dead load. As can be seen from the figure, the effect of dead loading is to improve the 
accuracy of the shapes at small loads. This is because, in the absence of any loading, 
the flexible sheet does not make a contact with all the indentors. The application of 
dead loading is to bridge the gaps left between the sheet and indentors, if any. With 
further increase in dead load the accuracy deteriorates. The effect of dead load on the 
RMS error is shown in Figure 3.10. 

Figures 3.11 and 3.12 present the variation of AVAB and RMS errors with an increasing 
UDL. The effect of UDL on the surface shape accuracy has been found to be more or 





RMS Error (mm) 
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Dead Load (kPa) 


Figure 3.10. Variation of RMS Error with Dead Loading. 
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less same as that of dead loading. This is evident from Figures 3.1 ] and 3.12. Based on 
the similar numerical experiments carried out with other surface shapes, the variation 
of AVAB and RMS errors with dead loading/UDL has been found to be qualitatively 
sanu as here. However, the load values up to which the shape accuracy improves, 
varied from shape to shape. 


3.5 Effect of Sheet Variables on the Surface Shape 
Accuracy 

The dimensions of the flexible sheet, the sheet thickness and the material of the flexible 
sheet form some of the important variables in the present study of shape realization. 
The effect of these parameters on the accuracy of realized shapes is important and 
helps in knowing the optimum values of these parameters for a given application. 

The dimensions of the sheet are important as it affects the realizability of a given shape 
as well as the size of the flexible surface tooling machine. The flexible sheet used here 
is of rectangular shape. The size of the sheet is larger than the area in which the 
indenlors are placed (Fig. 3.13). The extra area of the sheet which does not come in 
contact with the indentors is an important variable and can be measured in ter ms of 
margins. The effect of sheet margins on the surface shape accuracy for a typical case 
is given in Table 3.7. The indentors were placed in an area of 1.0 m by 1.0 m in this 
case. 


X and Y Margin 

AVAB Error 

RMS Error 

200.0 mm 

0.21750 

0.26821 

300.0 mm 

0.23716 

0.27479 

400.0 mm 

1.70882 

2.00238 

500.0 mm 

4.58881 

4.98410 


Table 3.7 Effect of Sheet Margins on the Surface Shape Accuracy. 

As can be inferred from the table above, better shape accuracies can be obtained at 
smaller margins of the flexible sheet. This is due to higher tensions of the flexible sheet 
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assoriiilrij with smalliT margin vahirs. At higher margin values the flexible sheet looses 
rontarl wit li murr and more indentors giving rise to poor realized shapes. 

Figures .’1. 1 1 and -1 1 5 show t he effert of sheet thickness on the surface shape accuracy in 
terms of A\’AB and RMS error.s respectively. In the present analysis the flexible sheet 
is idealized as a membrane. Because of this the the sheet thickness has no effect on the 
deformed configuration of the membrane when there is no external loading. However, 
the effert of shwt thickness in the presence of an external loading is significant. This 
is evident from the above figures which show the variation of shape error with sheet 
thickness in the presence of an UDL of 0.25 kPa. 

The variation of the error l>etween the desired and deformed surface shapes with the 
change in material properties of the flexible sheet are presented in Figure 3.16. The 
the ratio of first and second Mooney constants Ci/C^ is generally used to characterize 
nibber-like materials. The figure presents the variation of AVAB and RMS errors with 
increasing values of CxfCi. As the material of the flexible sheet tend to be more elastic 
with iin rea-siiig values of C\fC% better results were obtained at these values. 
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liidentor 



Figure 3.13. Sheet Margins. 
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Figure 3.14. Variation of AVAB Error with Sheet Thickness. 



RMS Error (mm) 
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Figure 3.15. Variation of RMS Error with Sheet Thickness. 
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ape Accuracy. 



Chapter 4 


SURFACE DISCRETIZATION 
FOR THE FLEXIBLE SURFACE 
TOOLING 


4.1 Introduction 


Realizability of a shape by flexible surface tool has been discussed in the previous chap- 
ter. The proposed tooling system is based on the concept of approximating a given 
continuous surface with a discrete one. The discrete surface here is in the form of 
finite number of closely spaced rigid tool elements known as indentors. The indentors 
are of square/circular/hexagonal cross-section with spherical ends. They occupy fixed 
positions in the form of a two dimensional layout, such that, their axes are perpen- 
dicular to the middle surface of the flexible sheet in its undeformed position. Each 
of the indentors can be moved along their axes independently and can be locked in 
their respective positions. When all the indentors are moved by appropriate heights 
the spherical ends of all the indentors together constitute a discrete approximation of 

a given continuous surface. 
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Given a layout of indcntors and the indentor geometry, the heights by which the in- 
dentors are to be raised to approximate a given continuous surface can be evaluated 
trivially. This is done by first localizing the desired continuous surface appropriately 
with respect to indcntors and then raising each of the indentors incrementally till it 
makes a contact with the desired surface. The open-loop approach to shape realiza- 
tion discussed in the previous chapter is based on deciding the indentor heights in this 
manner. 

In the proposed tooling, the number of indentors used for approximating a continuous 
surface is finite as well as limited. In this context, the geometry of the individual 
indentors, the indentor density and the indentor layout are important factors in de- 
signing an effective discrete surface. The study of these parameters with an objective 
to improve the accuracy of realizable shapes is the main subject of this chapter. The 
specific objectives of this study are as follows. 

• A geometric study of the discrete approximation of a continuous mold surface. 

• To study the change in effective mold shape due to presence of flexible sheet. 

• To study the change in effective mold shape due to external loads on the sheet. 


4.2 Indentor Density 

One of the important parameters related to discrete approximation of a continuous 
surface is the number of indentors used or the indentor density. Theoretically, a con- 
tinuous surface can be approximated by an infinite discrete surface tool elements of 
infinitesimal size. In practice, the number of surface tool elements or indentors used 
for approximating a continuous surface is finite. A small number of indentors used 
for approximation often yields poor results. When the number of indentors used for 
approximation is too large it is difficult to manipulate and control them individually. 
In this context, the effect of indentor density on the accuracy of realizable shapes is an 
important aspect. 

The concept of discrete approximation to a continuous surface is made use of, in an 
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impression toy for recording negative impressions of a complicated curved geometry 
such as palm of a human hand. The toy consists of a matrix of large number of pins 
arranged in an area. An individual pin can be pushed and locked in the longitudinal 
direction. The object whose impression is to be recorded is placed over the matrix 
of pins and is pressed against them gently until all the pins come in contact with the 
surface of the desired object. At this stage, the pins are locked in their respective 
positions and the object is removed. The end points of all the pins together form a 
discrete negative impression of the desired object. A discrete impression of a spherical 
object can be seen in Figure 4.1. 

In the case of a toy, the number of pins used is very large and the area of individual 
pins is very small. This can be considered as an example of approximating a continuous 
surface with discrete points. But, in the case of indentors the area of cross-section is 
significant and cannot be idealized as points. One way to measure the effectiveness of 
the discrete surface here is by the fraction of volume under the desired surface patch 
which is occupied by the indentors. For example. Figure 4.2 demonstrates this aspect 
for a two dimensional case. Here, a discrete approximation to a curve is built using 
area templates which are 2D equivalent of indentors. The figure shows the effectiveness 
of the approximation with increasing number of templates. 

This can be extended to a 3D case by considering a typical surface shape of compound 
curvature. Considering the surface patch given by the equation, 


z = 0.25 


1 - (5i/4)"] 


1 - (%/■*)“ 


where — 0.5 ^ ^ ^ 0.5 and 0.5 ^ j/ ^ 0-5 


(4.1) 


In order to approximate this surface patch, indentors in the range of 81(9 by 9) to 
225(15 by 15) are used. In each case the approximation error is plotted for number 
of indentors used (Fig. 4.3). The effectiveness of the discretization with increasing 
number of indentors is evident. The approximation error here is evaluated by the 

following expression. 
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Figure 4.2. A Discrete Approximation of a Curve by Area Templates. 



Approximation Error (%) 
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Number of Indentors 


Figure 4.3. Effectiveness of the Discretization with Increasing Number of 
Indentors. 
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. ■ , volume occupied by the indentors 

approximation error - 1 1- ^ ^4 2 ) 

volume under the surface patch 

In the proposed tooling, the discrete surface formed by a grid of indentors is not directly 
used as a mold surface. A flexible sheet is draped over the discrete surface to obtain a 
continuous surface. The presence of the flexible sheet over the indentors changes the 
effective mold shape. A physically based analysis accounting for the deformation of 
the flexible sheet is necessary to know the actual mold shape. Moreover, the presence 
of external loading on the flexible sheet further changes this shape. It is necessary to 
know the change in mold shape due to presence of external loading. 

Based on the theoretical formulation and the numerical scheme discussed in the sec- 
ond chapter, shape realization experiments were carried out with increasing number 
of indentors. The detailed procedure of shape realization has been discussed in the 
previous chapter. The results of one such experiment for the surface patch discussed 
above are given in Figures 4.4 and 4.5. in terms of AVAB and RMS errors respectively. 

In the absence of any external loading on the sheet, the AVAB and RMS errors contin- 
uously decrease with the increasing number of indentors. In the presence of an UDL 
on the sheet, there is a significant improvement in the accuracy of realized shapes in 
the range of 81(9 by 9) to 196(14 by 14). When the number of indentors are beyond 
this value the improvement in the accuracy is not much. It was pointed out in the 
previous chapter that with the application of external loading on the flexible sheet, the 
accuracy of realized shape first improves and then continuously decreases (Fig. 3.11 & 
3.12). It has been found that with the increasing number of indentors this tendency of 
the realized shapes to improve at smaller loads, decreases. 


4.3 Contact Geometry 

Contact geometry of the indentors is another important parameter in approximating a 
continuous surface for the tooling requirements. It refers to spherical end portion of the 
indentors which comes in contact with the flexible sheet. The contact geometry can be 



AVAB Error (mm) 
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Figure 4.4. Variation of AVAB Error with Increasing Number of Indentors. 




RMS Error (mm) 
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Figure 4.5. Variation of RMS Error with Increasing Number of Indentors. 
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measured m terms of radius of this spherical segment or contact curvature which refers 
to reciprocal of the contact radius. If the contact geometry of the indentor matches 
with the local shape of the surface patch to be realized, a better approximation can be 
expected. This is difficult to achieve assuming that all the indentors are of identical 
shape. Moreover, the same set of indentors are used for realizing a class of shapes 
varying significantly from one another. 

Considering the surface patch of previous section given by Equation (4.1), the effect 
of contact geometry on the realizability of a shape can be demonstrated. In the given 
parameter range, most of the region of this patch is synclastic. The variation of approx- 
imation error for this surface patch with increasing contact radii is shown in Figure 
4.6. The approximation error is evaluated using the expression given by Equation 
(4.2). The results of numerical experiments in terms of AVAB and RMS errors for this 
surface patch are given in Figures 4.7 and 4.8 respectively. For a given surface shape 
there exist a range of contact radius/curvature values of the indentors which gives best 
results. In this range the effect of external loading on the accuracy of realized shape 
is not much. In the absence of any external loading, best results were obtained with 
indentors of small curvature. This can be expected as the desired shape too has small 
values of curvatures. 

Similar numerical experiments were performed with other surface shapes of varying 
curvature values. Based on the results of these experiment it has been found that the 
same set of indentors cannot give good results for a class of shapes varying considerably 
in terms of curvature values. In order to improve the range as well as the accuracy 
of realizable shapes the use of removable indentors having different contact radii is 
suggested (Fig. 4.9). In this case, the indentors of required contact geometry can be 
inserted depending on the local shape of the surface to be realized. 


4.4 Indentor Layout 

In the proposed tooling, the arrangement of indentors or the indentor layout for shape 
realization forms another interesting topic. The indentor layout is closely related to the 
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radius = 4.9 radius = 6.53 radius = 8.17 



radius = 9.8 radius = 12.25 radius = 14.7 


pjgnj*g 4.9. Ins6rtabl6 Indantors of Diffarant Contact Radii. 
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indontor section geometry. Figure 4.10 shows four different combinations of indentor 
layout and indentor section geometry which are considered as a part of the present 
study. Due to compact arrangement of indentors, the layouts ’B’ and ’D’ with square 
and hexagonal sections respectively can be expected to give better results than layouts 
’A’ and ’C’. 

Table 4.1 presents the results of shape realization using four layouts for a typical 
shape. The contact geometry of the indentors in all the four cases is same. Based on 
the numerical experiments which were carried out on some typical shapes of varying 
curvature, the performance of layouts ”B” and ”D” was found to be better than layout 
”C”. The layout ”C” performed better than ”A”. The relative performance of the four 
layouts in the presence of uniformly distributed load of 0.25 kP a was also more or less 
same. Further, it is noticed that the variation in their performance was more in the 
case of realizing surface patches of small curvature values. 


Indentor 

Layout 

UDL 

(kPa) 

Average Absolute 
Error 

Root Mean Square 
Error 

Layout ’’A” 

0.0 

0.25 

0.9542 mm 
1.1665 mm 

1.2977 mm 
1.3132 mm 

Layout ”B” 

0.0 

0.25 

0.7978 mm 
1.0657 mm 

1.0551 mm 
1.2047 mm 

Layout ”C” 

0.0 

0.25 

0.8771 mm 
1.0869 mm 

1.0826 mm 
1.2351 mm 

Layout ”D” 

0.0 

0.25 

0.7429 mm 
0.9887 mm 

0.9500 mm 
1.1354 mm 


Table 4.1. Effect of Indentor Layout on the Surface Shape Accuracy 
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Figure 4*10. Indentor Layout. 





Chapter 5 


STATIC SHAPE OPTIMIZATION 


5.1 Introduction 

The results of shape realization experiments using flexible surface tooling were pre- 
sented in Chapters 3 and 4. In these numerical experiments, the relative positions 
of the indentors in order to realize a given surface shape were decided off-line. The 
procedure adopted for deciding the indentor heights has been discussed in the previous 
chapter. This is an open-loop approach to shape realization in which the error between 
the desired and realized surface shapes is fixed for a given set of tooling parameters. 
An alternate approach is to manipulate the heights of the indentors after every de- 
formation analysis run for any possible improvement in the accuracy of the realized 
shapes. The main objective of this chapter is to investigate this aspect in detail. 

The problem of shape realization by closed-loop approach can be formulated as a static 
shape optimization problem. This involves minimizing the error between the desired 
and realized surface shapes by treating the indentor heights as design variables. If the 
indentors are manipulated one at a time, the possible improvement in the accuracy 
of realizable shapes will be at the expense of high computational cost. Depending on 
the shape to be realized and the number of indentors used, the number of function 
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ovahiations may be excessively large. Such a scheme is not practical from the opera- 
tional i)oint of view and is not advisable in most of the cases. However, the number of 
function evaluations can be reduced by manipulating a group of indentors instead of 
manij)ulating a single indentor at a time. This aspect has been discussed in detail in 
the subsequent sections. 


5.2 Previous Work 

The problem of optimizing the static shape of a deformable object is often found in the 
area of design and control of large flexible space structures, such as, large antennas, 
mirrors, telescopes and interferometers. The shape control of these space structures 
is of great interest to structural designers (Haftka and Adelman, 1985; Orozco and 
Ghattas, 1991; Padula et al., 1989). These structures are often mounted on a number 
of discrete supports, some or all of which are fixed in space and are known as passive 
supports. If a support can be manipulated to change the static shape of the structure 
after erection, it is referred to as an active support. One of the main objectives here 
is to seek ways to minimize the need for active supports or to eliminate them by 
careful design and construction of the space structure. The studies in this area have 
shown that in some cases geometrical shape accuracy can be maintained with passive 
supports by carrying out off-line shape optimization. One of the ways to achieve the 
desired static shape of the reflecting surface here is to manipulate the supports without 
changing their respective locations. This is a case of finding the optimum supports for 
the structure when its topology and shape are predetermined. This problem is similar 
to the present closed-loop approach to shape realization by flexible surface tooling. 

The above problem is usually posed as a nonlinear optimization problem with an ob- 
jective to minimize the error between the desired and deflected shape of the structure. 
Conventionally, a structural analysis is first performed for some mitral guess of the 
design variables and using this as a starting point search direction is determined with 
the aid of a nonlinear optimization algorithm. The search direction is now used to 
improve the current design and the cycle is repeated until some convergence criteria is 
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met. This approach is usually called nested and several procedures have been proposed 
to improve its performance (Haftka, 1989; Orozco and Ghattas, 1991). 


5.3 Present Methodology 

A similar procedure is possible in the present case of static shape optimization of 
the flexible sheet for optimum indentor heights. The error between the desired and 
realized surface shapes can be minimized by using a nonlinear optimization technique 
and treating the indentor heights as design variables. The evaluation of objective 
function in this case involves the deformation analysis of the flexible elastomeric sheet 
for a given set of indentor heights. In the worst case, when there is no symmetry in the 
arrangement of indentors, the number of design variables here is equal to the number 
of indentors. Solving such a large scale optimization problem is highly prohibitive as 
it involves the evaluation of the objective function a few thousands of times by any of 
the first or second order optimization techniques. 

An alternate procedure which can be attempted in this case is to identify and manip- 
ulate a group of indentors which are likely to improve the accuracy of realized shapes. 
Based on the error between the desired and realized surface shapes at any particular 
time, a single or a group of indentors can be identified. This involves scanning the error 
profile and identifying the high error zones and the indentors in the vicinity of these 
zones. Once the indentors are identified their heights can be manipulated in propor- 
tion to the error in their respective zones. This is a case of an interactive optimization. 
Such an approach will involve only a few objective function evaluations. This cycle can 
be repeated until the error is more or less uniform throughout the flexible sheet. Once 
the indentors to be manipulated are identified the optimization in this case reduces to 
that of a line search. 

The accuracy of the realized shape can also be improved in some cases by moving all 
the indentors in either direction and keeping their relative configurations unchanged. 
Raising and lowering of all the indentors by an equal amount is equivalent to changing 
of the tension in the flexible sheet. Figure 5.1 shows this aspect for a two dimensional 
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(as( w.1k r( the flexible sheet is made to stretch continuously as the indentors are raised 
without changing their relative heights. 


5.4 Numerical Examples 

Three numerical examples are presented in this section. The first example deals with 
the shape realization of a doubly curved surface patch. The parametric representation 
of this surface patch is given as, 

^ = 0.15 [l-(5x/4)^] [l-(5y/4)^] (5.1) 

where —0.5 < x < 0.5 and —0.5 < y < 0.5 

In this case, first the deformation analysis of the flexible sheet has been carried out by 
following the open-loop approach discussed in the last chapter. A 4.0 mm thick sheet 
and 81 indentors are used which are arranged in the form of a matrix of 9 by 9. A 
uniformly distributed load of 0.5 kPa has been applied on the flexible sheet. Figure 
5.2 show^s the error profile between the desired and realized surfaces patches for this 
case. It can be seen from the figure that the error in this case is equally distributed 
throughout the flexible sheet. This indicates that the improvement in the realizability 
of the shape is not likely to happen by varying the relative heights of the indentors. 
In this case all the indentors can be raised and lowered by the same height to see the 
effect on the realizability of this shape. 

Figure 5.3 shows the effect of moving all the indentors on the accuracy of the realized 
shape in this case. It can be seen from Figure 5.3 that there exist an optimum height 
at which the best results can be obtained. At lower indentation heights, the tension 
in the flexible sheet is relatively small. In such case the effect of external loading on 
the realized shape is significant and responsible for relatively poor results. When the 
sheet is subjected to high tension values its ability to take the external load increases. 
However, at higher tensions, the tendency of the elastomeric sheet to form local areas 
of peaks and valleys increases (Fig 5.3). 
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Figure 5.2. Error Profile (Numerical Example I). 
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Th<- sc<o,id example is cenceriied with the realizability of another surface patch of 
ron.po„„<i curvature. The equation of the surface patch in this case is given as, 

a = 0.1 [(4a/3)^-l]"[(4j,/3f_l]“ (5.2) 

where -0.5 < a: < 0.5 and -0.5 < y < 0.5 

The d('formation analysis in this case is carried out with a total of 81 indentors and 
using a .5.0 nirn thick elastomeric sheet. An UDL of 0.5 kPa has been applied on the 
flexible sheet in an area where the indentors are placed. Figure 5.4 shows the error 
profile between the desired and realized shapes. The error distribution in the flexible 
sheet here is not uniform. This is a case of an example where the relative movements of 
the indentors may lead to the improvement in the shape accuracy. A group of indentors 
is identified in the region of high error values and the sensitivity of the movement of 
these indentors is first carried out. Based on this, a search direction is evaluated and 
the optimization is carried out using a quadratic search. Table 5.1 presents the results 
for this surface patch after two iterations of interactive optimization. The improvement 
in the accuracy of the realized shape is evident from results. 


Case 

Average Absolute 
Error 

Root Mean Square 
Error 

Open-Loop Approach 

1.2103 mm 

1.5077 mm 

Indentation Height Optimization 

1.0922 mm 

1.3949 mm 

After First Iteration 

0.8147 mm 

1.0267 mm 

After Second Iteration 

0.7931 mm 

1.0114 mm 


Table 5.1 Results of Static Shape Optimization (Numerical Example 11). 

The last example of this section is concerned with the realization of a swept patch 
given by the following equation. 

X = 1.51625 cos{u) 
y = l.O + u — 0.5 
2 = 0.758125 sm(u) - 0.658125 


(5.3) 




Chapter 5. Static Shape Optimization 


101 



Figure 5.4. Error 


Profile (Numerical Example II). 
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where 1.0534 < u < 2.0882 and 0.0 < t; < 1.0 

I he (ieforniatioii analysis is first carried out with a total of 81 indentors forming a 
mat nx of 9 hy 9. The layout of indentors in this case is shown in Figure 5.5. Figure 5.6 
shows the realized surface patch and the error profile between the desired and realized 
surface patches when there is an external load of 0.25 kPa in the form of UDL. 

1 he ( urvature of this surface patch varies only in one direction; the curvature in the 
direct ion of sweep being zero. All the indentors placed along the direction of sweep are 
raised by same height. This reduces the total number of design variables in this case to 
9. 1 he number of design variables further reduces to a value 5 as there is a symmetry 
in the shape to be realized and also in the placement of indentors. As the number of 
design variables here is small, the static shape optimization can be carried out in an 
automated manner. The conjugate gradient method(with Polak-Ribiere correction) is 
used for this optimization. The results of the same after two iterations are presented 
in Table 5.2. 


Case 

Average Absolute 

Root Mean Square 


Error 

Error 

Indentation Height Optimization 

0.7990 mm 

1.1024 mm 

After First Iteration 

0.7549 mm 

1.0606 mm 

After Second Iteration 

0.7056 mm 

0.9511 mm, 


Table 5.2 Results of Static Shape Optimization (Numerical Example III). 

The relative performance of the three schemes of static shape optimization discussed 
above is found to vary from shape to shape. In the first two examples considered above, 
the RMS error between the desired and realized shapes have been minimized whereas, 
the average absolute error forms the objective function in the third case. In all the 
three cases, a relatively small number of indentors are used. This is done intentionally 
to reduce the excessive CPU time involved in these computations and also with an 
objective to reduce the number of indentors required for the realizability of the typical 
shapes discussed above. 
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Realized Shape 



Error Profile 


Figure 5.6. Realization of a Swept 


Surface Patch (Numerical Example III). 


Chapter 6 


FLEXIBLE SURFACE TOOLING 
WITH BONDED INDENTORS 


6.1 Realizing Anticlastic Shapes 

A flexible surface tooling in which there is no bonding between the indentors and the 
elastomeric sheet was the subject of last three chapters. The realizability of a certain 
class of shapes has been demonstrated with this version of tooling. The tooling is found 
to be effective for synclastic shapes with small curvature values. This version of tooling 
can also be used for anticlastic shapes with very small Gaussian curvature values. For 
example, the numerical examples 3 and 5 of Chapter 3 dealt with shapes which have 
both regions of elliptic and hyperbolic points. However, the nature of tooling does not 
allow the shapes with considerable negative Gaussian curvature to be realized. Figure 
6.1 shows some typical shapes which cannot be realized with this version of tooling. 

In order to extend the applicability of the flexible surface tooling for shapes with 
considerable negative Gaussian curvature, a second version of the tooling has been 
proposed. Figure 6.2 shows a conceptual sketch of this version of tooling. It consists of 
a matrix of indentors whose ends are bonded to a flexible sheet. In the base position. 
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who,, » I the .ndenlors are at same height, the flexible sheet is in the undeformed 
state W hen the ,„dcntors are raised by different heights, the sheet deforms to take a 
speo,f„- shape. The sheet can be made to deform to a desired shape by manipulating 
t he mdentor he,ghts. A numerical study of the shape realizability by this version of 
tooling IS the main subject of this chapter. 

Umlike tlie earlier version of tooling, the indentors-in this case are placed at some 
distance apart. This is necessary to allow the the portion of a sheet which is not in 
contact with the indentors to deform freely. As the indentors are placed a distance 
apart, their slenderness ratio is kept as low as possible to avoid any possible buckling 
due to reactive contact forces. The advantage of flexible surface tooling with bonded 
indentors is that it can be used for realizing shapes having positive or negative Gaussian 
curvature or both. However, the nature of tooling restricts the realizable shapes to 
those having small curvature values. This is due to flexible sheet being subjected to a 
high local straining in the regions of large curvatures. 

The indentors can be bonded to flexible sheet either by some mechanical means or 
using suction jacks. When suction cups are used the process is similar to that of 
thermoforming. With a mechanical bonding, rotational degrees of freedom at the point 
of bonding can be incorporated. For the purpose of analysis, the contact between the 
the indentors and the flexible sheet can be treated as a point contact. The deformation 
analysis of the flexible sheet in this case involves finding of equilibrium shape of the 
sheet for certain positional constraints. The shape realization program described in 
Chapter 3 has been modified to account for this new boundary condition. 

6.2 Numerical Examples 

Five numerical examples are presented in this section. In order to compare this version 
of tooling with the earlier one, first two examples here are taken from that of Chapter 
3. The first example is a linearly swept surface patch whose generatrix is an elliptic 
arc and the directrix is a line normal to the plane of generatrix curve. The parametric 
equation of the patch is reproduced here for the sake of convenience. It is given by. 
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^ — 1.51625 cos(u) 
y = t;-0.5 

2 = 0.758125 6m(u) - 0.658125 (6.1) 

where 1.0534 < u < 2.0882 and 0.0 < n < 1.0 

In Older to realize this patch an elastomeric sheet of 1.1 m by 1.1 m and thickness 
o.O 77int is used. The deformation analysis is carried out with a total of 121 indentors 
bonded to a flexible sheet. The layout of the indentors can be seen in Figure 6.3. The 
ratio of first and second Mooney constants, C\IC 2 is taken as 4.0. The analysis has 
been carried out with a total of 1936 triangular elements (Fig. 6.4). The results of the 
analysis in terms of AVAB and RMS errors are given in Table 6.1. Figure 6.5 shows 
the realized shape and the error between the desired and realized shapes at higher 
magnification. 


Indentor-Sheet 

Average Absolute 

Root Mean Square 

Interaction 

Error 

Error 

Bonded 

0.4817 mm 

0.5521 mm 

Unbonded 

0.6936 mm 

0.8611 mm 


Table 6.1 AVAB and RMS Errors (Numerical Example I). 


The second example is concerned with the realization of a doubly curved surface patch. 
The patch is defined by the following expression. 


z = 0.25 


1 — (5a: /4)^ 


1 - (5y/4)^ 


where —0.5 ^ x ^ 0.5 and —0.5 ^ y ^ 0.5 


( 6 . 2 ) 


The analysis in this case is performed with a 4.0 Tfitn thick sheet. The other input data 
for this analysis is same as that of previous example. The error between the desired 
and the deformed patches in this case is given in Table 6.2. 
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Figure 6.3. Indentor Layout (Numerical Example I). 




Bonded rnH-ntn,. 
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Figure 6.4, Finite Element Discretization (Numerical Example I). 


patch 
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Indentor-Sheet 

Interaction 

Average Absolute Root Mean Square 
Error Error 

Bonded 

0.9911 mm 

1.3243 mm 

Unbonded 

1.1948 mm 

1.4755 mm 


Table 6.2 AVAB and RMS Errors (Numerical Example II). 

For the two surface patches discussed above better results were obtained with indentors 
I fonded to flexible sheet. However, the results in practice may diflFer slightly from the 
reported values due to contact between the flexible sheet and the indentors being not 
a point contact. 

The next two examples deal with the realization of shapes which have both synclastic 
and anticlastic regions. In the third example the results for a linearly swept surface 
patch are presented. The generatrix of the swept patch is a sinusoidal curve. The 
parametric equation of the patch is given as, 

X ~ u 
y = V 

z = A sin{2'!ru) (6-3) 

where 0.0 ^ u ^ 1.0 and 0.0 ^ v ^ 1.0 

In this case, the analysis has been carried out for various values of A in the range of 
0.02 to 0.2 m using a 5.0 mm thick sheet with 169 indentors. The realized surface patch 
and the error between desired and realized patches are shown in Figure 6.6. Figure 6.7 
shows the variation of AVAB and RMS errors for this patch for various values of A. At 
smaller values of A which corresponds to surface patches with small curvature values, 

better results were obtained. 

The next example is concerned with the realisation of a bicubic Bezier patch. The 

parametric equation of the patch is given as, 

X = n — 0.5 
y = V — 0.5 




4 




Error (mm) 
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Figure 6.7. Variation of Realized Shape Accuracy with Amplitude ”A 

(Numerical Example III). 
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2 = 20uu(l -u)(l-i;)[-10 + 30u + 25i;-70H (6.4) 

where 0.0 < u < 1.0 and 0.0 < v < 1.0 

The results of the analysis for this patch are given in Table 6.3 and the realized shape 
is depicted in Figure 6.8. The error between the desired and realized surface patches 
at higher magnification is also shown in Figure 6.8. It can be seen from this figure that 
the distribution of shape error is not uniform throughout the flexible sheet. The error 
is more in the regions of high curvature. 


Average Absolute 

Root Mean Square 

Error 

Error 

0.8861 mm 

1.2718 mm 


Table 6.3 AVAB and RMS Errors (Numerical Example IV). 


The last example of this section is a hyperbolic paraboloid patch. This is an example 
of anticlastic surface patch. The parametric form of this surface patch is given below. 


p(u,t)) = [1 - U U 


0.0 0.0 0.00 

1.0 0.0 0.24 


T 

T 


0.0 1.0 0.24 

1.0 1.0 0.00 


T " 
T 



(6.5) 


where 0.0 ^ ^ 1*0 ^nd 0.0 ^ ^ ^ 1*0 

The results for this surface patch are given in Table 6.4. The realized patch is shown 
in Figure 6.9. The results obtained for this patch are much superior than some of the 
earlier cases discussed above. Further numerical experiments related to realization of 
hyperbolic paraboloid patches which were carried out produced similar results. The 
better results obtained in the case of hyperbolic paraboloid patches can be explained 

in the following way. 





Error 


profit® 



Figure 6.9. 
ample V). 
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Average Absolute 

Root Mean Square 

Error 

Error 

0.0486 mm 

0.1158 mm 


Table 6.4 AVAB and RMS Errors (Numerical Example V). 

The deformed configuration of the flexible sheet in this version of tooling mainly de- 
pends on the displacement constraints imposed by the indentors. Considering a set of 
four adjacent indentors in the form of a rectangular grid, a boundary can be drawn by 
joining the four end points of these indentors by straight lines. An infinite number of 
surface patches can be constructed interpolating these four boundary edges. Among 
all these surface patches, the one which has a minimum surface area is the hyperbolic 
paraboloid defined by four end points of the indentors. As the rubber-like membrane 
has a tendency of deforming in a form minimizing its deformed surface area, supe- 
rior results were obtained in the case of hyperbolic paraboloid (Callahan et al, 1988; 
Grundig, 1988; Hansen, 1993). 

In the numerical examples discussed above, it is assumed that there is no external 
loading on the flexible sheet. The numerical experiments were carried out to study the 
effect of uniformly distributed load (UDL) on the surface shape error. The effect of 
external loading on a typical realized shape is shown in Figure 6.10. The error between 
the desired and realized surface shapes continuously increases with an increasing loa . 


6.3 Study of Tooling Parameters 

The effect of tooUng pa^enteters o„ the accuracy of reatable shapes h^ b^n discussed 
in Chapters 3 aud 4 for the flexible surface tooling in »hich there is no bonding between 
the flexible sheet and the indentors. Not aU the tooling 
relevant in the present c^. The effect 

of readable surface shapes are presented in this section, tne p 
here are the sheet thichness, the sh,*t material and the indentor density. 




Surface Shape Error (mm) 
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Figure 6.10. Effect of UDL on the Surface Shape Accuracy. 
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Tlu' shape realization experiments were performed with different sheet thickness in the 
range of 2.0 mm to 10.0 mm. Figures 6.11 and 6.12 present the effect of sheet thickness 
on tlu' surface shape accuracy. In the absence of any external loading, the thickness 
does not effect the shape accuracy. This is due to membrane state of the flexible sheet. 
When there is an external load in the form of UDL acting on the sheet, the accuracy 
of realized shapes continuously increases with the increasing sheet thickness. The use 
of higher thickness of the sheet is advisable when there is a considerable external force 
acting on the flexible sheet. However, as the thickness increases the the reactive contact 
forces on the indentors bonded to flexible sheet also increases. A component of this 
recative contact force acting transversely on the indentor may lead to possible buckling 
if the slenderness ratio of the indentor is not high. 

The effect of sheet material on the realizability of a typical shape is showm in Figure 
6.13. The figure shows the variation of AVAB and RMS errors for an increasing values 
of C-ifC 2 , where Ci and C 2 are the first and second Mooney constants respectively. 
Unlike the earlier version of a flexible surface tooling discussed in Chapter 3, the effect 
of sheet material on the realizability of a shape is not much in the present case. 

The number of indentors used for shape realization is another important factor which 
has been considered for the present study. Indentors in the range of 81(9 by 9) to 
289(17 by 17) are used for the analysis. In each case the AVAB and RMS errors are 
plotted for number of indentors used (Fig. 6.14). The accuracy of realized shapes 
strongly depends on the number of indentors used. This is also true with the earlier 
version of the tooling discussed in Chapter 3. 
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Figure 6.12. Effect of Sheet Thickness on RMS Error. 
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Chapter 7 


CONCLUSIONS 


7.1 Summary 

In th(' present work, a flexible surface tooling system has been proposed for realizing 
free-form shapes. The tooling is based on the concept of discrete approximation to a 
cont.inuous mold surface. This is a variation of the discrete surface tooling proposed 
earlier mainly for sheet metal forming. The applicability of the proposed tooling has 
been st udied for shape conforming processes where a single finished surface is desired 
such as composite layup. Unlike the earlier investigations in this direction, which were 
based on laboratory experiments, the present approach is a computational one. 

A computer code for the deformation analysis of rubber-like membranes in multiple 
contact has been developed for carrying out the proposed theoretical study. The effect 
of various parameters related to the proposed tooling on the realizability of a shape has 
been studied in order to identify a typical range of these values for a given application. 
The realizability studies have been carried out by two different approaches namely 
open and closed-loop approaches which are based on deciding the indentor heights in 
an off-line and on-line manners respectively. The range of realizable shapes by the 
proposed tooling include both synclastic and anticlastic shapes. 
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• Tlu' accuracy of realized shape can be considerably improved in some Ccises by 
following a closed loop approach which involves the manipulation of a single, a group 
or all the ind('ntors after every deformation analysis run. 

• The flexible surface tool in which the indentors are bonded to flexible sheet can 
b(' usc'd for realizing anticlastic shapes with small curvature values. 


7.3 Scope for Further Work 

The present work is concerned with a numerical study of the shape realization using the 
])ropo.s('d flexible surface tooling system. The concept can be extended to all the shape 
realization processes with appropriate modifications. Though the results presented 
here are quantitative in nature, the design inferences are more qualitative in nature. 
It has been felt that the numerical study of the proposed flexible surface tool can be 
improved by incorporating one or more of the following features. 

1. Accounting for the contact friction between the flexible sheet and the indentors. 

2. A viscoelastic analysis of the deformation process by including the thermal effects 
which are likely to be present in the case of processes like Layup. 

3. Carrying out the deformation analysis with 3D finite elements or 2D elements 
with bending stiffness. 

4. Extending the analysis to calculate the stresses in the flexible sheet and the 
possible deflections of the indentors in case of tooling with bonded indentors. 

Most of the relevant tooling parameters which affect the realizability of a shape in the 
proposed tooling are covered as a part of the present work. However, a more detailed 
study of the these parameters using the concept of design 0 / experiments is necessary 
in this case. The need tor further numerical experiments to establish a close rela ion 
between the tooling parameters and the complexity of the shape to be realized was also 
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h'lt. I’lif' j)r()r('.ss of sliape optimization discussed in Chapter 5 which was carried out 
int('ra('t ivcly here, needs to be automated using an appropriate shape control algorithm. 

In tlu- pr<'sent work of shape realization, no attempt has been made to compare the 
n'sults of computational analysis with the physical experiments. This is necessary to 
get a more clear picture of the realizability of a shape by the proposed tooling. The 
effort .s to build a flexible surface tooling system are already in progress and some of the 
If'ssoiis h'arnt from the computational analysis are being incorporated in this endeavor. 


7.4 Scope for Related Work 


Some of the concepts and techniques which are the outcome of the present work can 
bo used in other areas of engineering. Two such applications which have similarities 
with the present problem of shape realization are geometric modeling of draped fabric 
surfaces and the shape finding of large spatial structures. 


Geometric modeling of fabric-like objects has been an active area of research for the 
past one decade. The fabric-like materials when draped on support geometries, deform 
int o complicated doubly curved surfaces which are aesthetically pleasing (Collier, 1990-, 
Dowlen, 1976). The objective here is to model the shape of the deformed fabric when it 
is draped on the complicated support geometry. In another version of the problem the 
objective is to achieve a close conformance between the the fabric and the rigid surface 
on which it is draped (Aono et ah, 1994-, Hinds et al, 1991). These problems have 
applications in apparel design, fashion technology, computer animation and computer 

synthesized advertisements. 


Initial attempts to model fabric-like objects were pure geometric in nature which . 
not yield realistic results. On the other hand, a finite element model accounting or t e 
complex deformation process of the fabriclike materials is computational y intensive 
in nature (Lloyd, 1984). Thus attempts have been made 

by hybrid geometric-physical deformable models (Coquillarl, 1990; ^e^ 

nopoulos, 1988; Thingvold and Cohen. 1990; Weil, 1986). The results of one such 
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model proposed by the author for fabrics is shown in Figure 7.1 (Dhande et ah, 1993). 

I 1h' problem of finding the shape of a draped fabric can also be posed as a mechanics 
pr<.d)!em similar to the present Ccise of deformation analysis of rubber-like flexible sheet. 
It can be further simplified by making some of the assumptions about the deformed 
geometry which are often used in the area of textile mechanics. Such an attempt has 
b(>en made to model fabric-like objects in parallel to the present research work (Rao 
and Dhande, 1995). 

Another area of research where some of the concepts proposed in the present work have 
aijplirations, is the shape finding and shape optimization of the large shells and spatial 
structures (Berger, 1988; Fujikake et ah, 1989; Ramm, 1992; Sclaich et ah, 1989). The 
two geometric forms which are often used for large spatial structures such as roofs of 
large auditoriums and stadia are fabric membrane structures and pneumatic membrane 
structures. The first step in the design of these two forms of tensile membrane struc- 
tures is the process of shape finding. Shape finding often involves deforming a flexible 
membrane under a controlled load case like pressure or dead load to a desired shape 
looking for a pure membrane stress state. It is common to use analytical, numerical 
and experimental methods for the process of shape finding. Often the shape of the 
spatial structure to be built is realized in a laboratory using models specifically built 
for this purpose (Aucouturier, 1966; Isler, 1993; Kirsh and Reiss, 1967). It is felt that 
a generic shape realization tool for realizing a variety of shapes is still lacking. Some 
of the lessons learnt from the shape realization experiments using flexible surface tool 
of the present work will be useful for building one such shape realization machine. 
Moreover, the present computational analysis too has many things in common with 
the shape finding of pneumatic membrane structures by numerical methods such as 

FEM. 



Figure 7.1. A Hybrid Geometric-Physical Deformable Model for Fabrics. 
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